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Introduction 

in  many  cases  vortices  are  being  generated  on  the  aircraft 
surface:  on  the  upper  wing  and  body  surfaces,  in  the  vrcinr- 
tv  of  the  deflected  fins.  They  can  essentially  influence 
pressure  distribution  and  heat  transfer.  The  vortex  genera¬ 
ted  on  the  upper  surface  of  a  body  or  a  wing  can  encounter 

the  shock  wave  in  front  of  the  fin,  located  downstream 

1 1  e  •  T pTSLct i on  of  th-B  vortsx,  9^ 
fFiq.la).  Another  example  is.  interaction  u  .  ,  ^  . 

^  a  cxtrake  or  a  nose  flap  intended 

nprated  by  a  canard  [1],  a  strake  or 
neratea  y 

for  aircraft  trimming  with  the  wing 

(Fig  lb).  A  case  important  for  practical  purpose  is  descri¬ 
bed  in  [2]  (Fig.lc):  Interaction  between  the  vortex  and  in 

let  shock  waves  qualitatively  changes  the  inlet  entry  flow 
pattern.  A  similar  phenomenon  can  be  observed  in  ° 

interaction  of  the  vortex  generated  in  the  inlet  wi 
pylon  of  the  air-breathing  jet  combustion  chamber  (Fig.l  ). 
Tn  this  context  it  is  worthwhile  to  Investigate  the  vortex 
structure  and  its  interaction  with  the  streamlined  body. 

Wing  flow  at  an  angle  of  attack,  followed  by  a  flow  separa¬ 
tion  and  vortices  formation,  have  been  investigated  in  many 
studies.  in  some  papers  vortices  structure  is  investigated 
in  detail.  Most  works  deale  with  flows  of  incompressib  e 

fluid  A  vortex  layer,  shedding  from  a  wing  surface,  ro  s 
up  into  a  spiral.  Its  shape  and  velocity  distribution  are 
approximately  described  by  .the  equations  of 

compressible  fluid  [3,4].  At  high  Reynolds  numbers  the  expe¬ 
rimental  [5,61  and  theoretical  results  are  in  9°°^ 
ment.  However  velocity  tangential  discontinuity  on  the  spi 
ral  surface  typical  for  an  incompressible  fluid  flow  ar 
smoothed  due  to  the  influence  of  viscous  di  us 
The  influence  of  viscousity  is  still  greater 

vortex  core.  It  was  discovered  by  the  experiment  [5]  and  n 
vestlgated  theoretically  for  a  laminar  flow  in  works  [7,8]. 
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A  survey  of  these  investigations  is  given  in  [9].  In  1957 
peokham  and  Atkinson  discovered  a  vortex  burst  [10).  They 
observed  it  in  the  flow  field  in  the  vicinity  of  the  wing  in 
the  region  of  a  positive  pressure  gradient.  The  vortex  burst 
is  also  observed  in  the  swirling  nozzles  and  diffusers 
flows.  After  the  discovery  of  the  vortex  burst  a  great  vari¬ 
ety  of  theoretical  investigations  of  this  phenomenon  was 
carried  out. 

A  supersonic  vortex  structure  was  thoroughly  investigated  in 
the  work  [1].  Using  modern  measurement  techniques  the  rese¬ 
archers  obtained  the  data  concerning  the  influence  of  the 
swirl  strength  and  pressure  gradient  on  the  velocity  and  gas 
pressure  in  the  vortex  at  Mach  numbers  from  1.6  to  2.28. 

In  1975  research  [2]  discovered  a  vortex  breakdown  in  front 
of  the  inlet  due  to  the  influence  of  a  strong  shock  wave. 
The  investigation  was  carried  out  at  Mach  number  from  1 . 4  to 
3.0.  It  was  shown  that  a  vortex  breakdown  is  followed  by 
formation  a  dead  region  which  induces  a  cone  shock.  The  inc¬ 
lination  angle  of  the  dead  region  borderline  and  pressure 
increase  in  it  are  close  to  similar  values  for  a  separation 
region  of  the  turbulent  boundary  layer.  The  formation  of  the 
dead  region  reduces  the  pressure  recovery  coefficient  in  the 
supersonic  inlet.  This  investigation  was  continued  in  works 

[1,11,12]. 

Paper  [1]  determines  the  paraiileters ,  influencing  the  vortex 
burst  when  interacting  with  the  normal  shock:  vortex  streng¬ 
th  and  Mach  number  (shock  pressure  difference).  The  vortex 
breakdown  conditions  are  specified  for  Mach  number  up  to 
2.5.  In  this  paper  on  the  basis  of  the  Euler  equations  are 
also  theoretically  determined  the  conditions  under  which  the 
vortex  breakdown  occurs.  The  calculation  results  and  the  ex¬ 
perimental  data  are  in  good  agreement.  In  paper  [11]  similar 
analysis  of  interaction  between  a  vortex  and  an  oblique 
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shock  past  which  .a  supersonic  flow  is  retained  was  carried 
out  at  Mach  numbers  of  2.28  and  5.0.  In  calculations  a  vor¬ 
tex  breakdown  was  not  obtained.  In  paper  [2]  it  is  shown 
that  at  Mach  numbers  from  2.5  to  4.0  a  vortex  breakdown  oc- 
ours  at  a  rather  high  angle  of  attack  of  the  wing  generating 
,  vortex  (at  a  >  5-10“).  In  paper  [12]  a  vortex  burst  was 
also  observed  during  Interaction  between  a  vortex  and  a 
central  shock  in  the  under-expanded  jet.  It  occurred  at  ang¬ 
les  of  attack  of  the  wing-vortex  generator  greater  than  8 

in  the  known  papers  a  stagnation  temperature  distribution  in 
a  supersonic  vortex  was  not  investigated.  At  the  same  time 
measurements  performed  at  subsonic  speeds  in  the  Ranque  tube 
Indicate  that  vortex  temperature  can  essentially  change 
[13] .  This  property  of  vortex  flows  is  applied  in  the  simp- 

lest  refrigerators. 

None  of  the  studies  familiar  to  the  authors  deals  with  the 
investigations  of  the  vortex  influence  on  the  heat  exchange 
on  the  frontal  surfaces  of  a  body.  At  the  same  time  the  ob¬ 
tained  data  about  the  vortex  breakdown  indicate  that  this 
influence  can  be  essential. 


The  aim  of  this  paper  is  as  follows: 

1)  To  obtain  additional  data  about  the  vortex  structure  in 
supersonic  flow  by  means  of  experiment  and  theoretical  in- 

vest igations . 

2)  To  obtain  experimental  data  about  the  vortex  influence  on 
the  heat  transfer  on  the  frontal  surface  of  the  blunted  bo¬ 
dy. 
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Nomenclature  for  the  part  A 


=  thermodiffusion  coefficient, 

=  chord  of  vortex  generator, 

=  heat  capacity, 

=  vortex  diameter , 

=  reference  vortex  diameter  (at  the  begin  of 
the  vortexes  straight  portion), 

=  sphere  diameter , 

=  at /6^ /  Fourier  number, 

=  q/(To-Tw).  heat  transfer  coefficient, 

=  enthalpy, 

=  length  of  separation  region  in  front  of  sphere, 

=  Mach  number, 

=  pressure, 

=  heat  flux, 

.  theoretical  value  of  the  heat  flux  in  the  stagnation 
point  of  the  sphere, 

=  D/2,  sphere  radius, 

=  u^  D  /  ,  Reynolds  number, 

=  Reynolds  number  (per  1  m)  , 

=  dTstrnce  along  sphere  surface  from  its  center, 

=  temperature, 

=  temperature, 

=  velocity, 

=  coordinate  in  free  stream  direction, 

=  distance  from  trailing  edge  of  vortex  generator 

to  sphere, 

=  coordinates  in  the  directions  perpendicular  to  a  free 
stream, 

=  angle  of  attack  of  a  vortex  generator, 

=  inclination  angle  of  shock  wave, 

=  azimuthal  angle, 

=  model  wall  thickness, 

=  angle  of  borderline  inclination  of  separation  region. 
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thermal  conductivity  coefficient, 
viscosity  coefficient, 
density, 
time, 

angle  between  the  x-axis  and  the  probe  axis 

Subscripts 

eo  =  free  stream  conditions, 

0  =  wind  tunnel  settling  chamber, 

i  =  initial  conditions, 

m  =  maximum  value, 

s  =  stagnation  point  past  normal  shock  wave, 

w  =  on  the  model  surface. 


I  = 
M  = 
P  = 

X 
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A.  Experimental  investigation  of  the  vortex  structure 
and  its  interaction  with  a  sphere 

1.  Experimental  facility 

1.1.  Description  of  the  T-33  wind  tunnel 

The  T-33  facility  is  an  intermittent,  blowdown  wind  tunnel 
(Fig. 2).  It  is  often  used  at  TsAGI  for  investigation  of  heat 
transfer  and  pressure  distribution  on  the  bodies  surface  in 
supersonic  gas  flows  as  well  as  for  flow  visualization.  This 
wind  tunnel  was  being  used  for  similar  investigations  in  the 
subsonic  flow. 

The  scheme  of  the  T-33  wind  tunnel  is  shown  in  Fig. 3.  In  the 
majority  of  experiments  an  electrical  resistance  heater  (1) 
is  used  for  air  heating.  Its  power  is  up  to  1000  kW.  When 
the  experiment  purpose  is  a  flow  visualizaton  another  heater 
is  used  wherein  heating  occurs  due  to  combustion  of  kerose¬ 
ne.  It  represents  a  modified  combustion  chamber  of  the  tur¬ 
bojet  engine. 

The  maximum  temperature  of  the  air  in  both  heaters  is  about 
500° C.  However  in  the  majority  of  experiments  the  temperatu 
re  doesnot  exceed  200*^ C.  In  settling  chamber  (2)  the  flow  is 
straightened  by  a  honeycomb  and  a  system  of  grids.  The  maxi¬ 
mum  pressure  in  the  settling  chamber  reaches  8*10®  Pa.  For 
developing  a  supersonic  flow  axisymmetric  nozzles  (3)  conto 
ured  to  produce  uniform  Mach  3,  4  or  5  flow  are  used.  There 

is  also  a  convergent  nozzle  generating  a  subsonic  flow.  In 
this  investigation  the  nozzle  intended  for  Mach  number  of  3 
was  used.  The  nozzle  exit  section  diameter  is  304.5  mm. 

The  wind  tunnel  test  section  (4)  (see  Figs.  3  and  4)  pre¬ 
sents  itself  the  Eiffel  chamber.  Its  diameter  is  800  mm,  the 
distance  between  the  nozzle  and  the  diffuser  is  550  mm.  In 
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the  walls  of  the  test  section  there  are  5  windows  for  illu¬ 
mination  and  observation  of  models.  For  fixing  the  models 
there  are  2  supporting  mechanisms  in  the  wind  tunnel  test 
section  (above  and  below).  They  can  inject  a  model  and  a  me¬ 
asuring  unit  or  two  parts  of  a  compound  model  into  the  flow 
simultaneously  or  one  after  the  other.  The  lower  supporting 
device  can  be  moved  in  the  axial  direction  during  the  expe¬ 
riment.  Furthermore  it  is  possible  to  change  the  model  angle 
of  attack  during  the  experiment. 

Below  on  the  test  section  wall  a  thermostatic  chamber  is  mo¬ 
unted.  During  the  heat  transfer  investigations  this  chamber 
prevents  the  model  installed  on  the  lower  supporting  mecha¬ 
nism  from  heating  at  the  period  when  the  flow  becomes  stea¬ 
dy.  After  obtaining  a  prescribed  steady  regime  a  flap  of  the 
thermostatic  chamber  is  moved  aside  and  the  model  is  quickly 
injected  into  the  flow.  A  maximum  model  length  is  350  mm. 

For  obtaining  the  required  air  rarefaction  in  the  test  sec¬ 
tion  a  two-stage  supersonic  ejector  is  applied. 

A  great  number  of  investigations  of  heat  transfer  by  means 
of  thermosensitive  coatings  and  discrete  thermocouples  was 
carried  out  in  the  T-33  wind  tunnel.  Furthermore,  a  flow  vi¬ 
sualization  by  a  shadowgraph  technique,  a  laser  knife  and  a 
calibrated  oil  points  is  often  performed  in  the  tunnel.  Some 
results  obtained  in  the  T-33  wind  tunnel  are  presented  in 
papers  [14,15]. 


1.2.  Vortex  generators 

Different  ways  of  vortex  generation  are  known.  In  a  number 
of  works  for  generation  of  a  vortex  a  swept  or  a  strainght 
wing  at  an  angle  of  attack  is  used.  Delery  [1]  used  for  vor¬ 
tex  generation  in  subsonic  flow  a  wing-type  generator  along 
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with  a  vane-type  generator.  The  vane-type  generator  presen¬ 
ted  itself  an  axi symmetric  body  on  which  two  vanes  were  ins¬ 
talled  at  a  certain  angle  of  attack.  In  Ref. [16],  a  vortex 
produced  also  by  swirl-vanes  being  placed  into  axisymmetric 
casing  ,  is  described.  The  investigation  [1]  has  shown  that 
the  velocities  distributions  in  the  vortices  produced  by  a 
wing  and  swirl-vanes  differ  substantially.  In  a  supersonic 
flow  vanes  without  external  casing  evidently  were  not  used 
up  till  now. 

In  the  current  investigation  three  generators  of  both  type 
were  used:  2  generators  in  the  form  of  a  rectangular  wing 
and  one  generator  of  a  vane  type.  In  designing  the  wing-type 
generator  the  investigation  results  of  Ref. [17]  were  used. 
In  this  paper  information  about  vortex  dimensions  and  its 
position  at  different  wing  angles  of  attack  was  obtained. 
Figures  5,  6  and  7  show  the  schemes  of  the  generators  which 

were  used.  Generators  1  and  2  have  an  rectangular  shape.  The 
leading  and  lower  edges  are  sharp  and  the  trailing  edge  is 
blunt.  Over  a  greater  part  of  the  wing  span  a  cross  section 
profile  has  a  triangular  shape.  Generators  1  and  2  differ 
only  by  the  chord  length.  The  wing  chord  of  generator  1  is 
almost  3  times  greater  than  that  of  generator  2.  Generator  1 
is  designed  mainly  for  investigating  a  total  pressure  and 
total  temperature  distribution  in  a  vortex  and  generator  2 
mainly  for  investigating  a  vortex  development  and  its  inte¬ 
raction  with  a  sphere. 

Generator  3  (Fig. 7)  presents  itself  an  axisymmetric  body 
with  a  diameter  of  6  mm,  on  which  4  rectangular  vanes  are 
fixed  at  an  angle  of  attack  of  30°  (there  are  also  genera¬ 
tors  with  the  vane  setting  angles  of  10°  and  20  but  they 
were  not  used) .  A  potential  advantage  of  generator  3  as  com 
pared  with  two  other  generators  consists  in  a  less  flow  per¬ 
turbance  because  the  plate  on  which  the  generator  is  mounted 
is  installed  at  a  zero  angle  of  attack. 
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1.3.  Total  pressure  and  stagnation  temperature  probes 
and  device  for  their  positioning 

A  stagnation  pressure  probe  (Fig. 8a)  presents  itself  a  ca 
pillary  tube  fabricated  from  a  stainless  steel  with  an  ex¬ 
ternal  and  internal  diameters  0.8  and  0.48  mm  respectively. 

A  capillary  with  a  straight  section  length  12  mm  was  bent  in 
a  manner  as  shown  in  Fig. 8a  and  went  through  a  positioning 
device  rod.  The  end  of  the  capillary  was  connected  with  a 
miniature  pressure  probe  with  a  measuring  range  from  0  to  10 
bar,  by  means  of  a  flexible  pipe.  The  inertia  of  the  total 
pressure  probe  did  not  exceed  0.1  s. 

A  stagnation  temperature  probe  (Fig. 8b)  presented  itself  a 
tube  made  of  a  stainless  steel  with  external  and  internal 
diameters  2  and  1 . 5  mm  respectively  and  straight  section 
length  20  mm.  Inside  this  tube  a  thermocouple  junction  from 
chromel-copel  alloys  with  a  wire  thickness  of  0.1  mm  was 
placed.  The  rear  end  of  the  tube  had  a  plug  and  there  were  2 
holes  0.5  mm  in  diameter  in  the  side  surface  behind  the 
thermocouple  junction.  A  signal  from  the  thermocouple  was 
recorded  by  a  data  acquisition  processing  system.  A  tempera¬ 
ture  recovery  coefficient  of  this  probe  design  at  specified 
flow  paramemeters  according  to  literature  data  was  estimated 

T 

—  0 . 965  ±  0 . 025 

T 

The  real  values  of  the  recovery  coefficient  were  estimated 
in  the  tests. 

For  moving  the  probes  a  positioning  device  was  developed  and 
fabricated  (Figs. 9, 10).  Probe  1  is  fixed  inside  the  rod  2 
(Fig. 9).  Miniature  cooled  engines  4  and  10  move  the  rod  in 
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the  vertical  <Iirection  and  rotate  it  relative  to  the  axis. 
The  positioning  device  is  mounted  on  the  lower  supporting 
mechanism  of  the  T-33  tunnel  and  can  be  moved  along  the  tun¬ 
nel  axis. 

The  positioning  device  allows  the  probe  to  be  moved  along 
the  y-axis  (in  vertical  direction)  at  the  distance  up-to  120 
mm  with  the  minimum  spacing  of  0.05  mm.  At  every  specified 
station  along  the  y-axis  the  positioning  device  can  stop  the 
probe  and  rotate  it  up  to  the  angle  ±90°  with  the  minimum 
setting  0.5°. 

The  positioning  device  is  controlled  by  a  computer .  A  cont¬ 
rol  code  allows  to  change  movement  and  rotation  steps  Ay  and 
A(p,  respectively,  as  well  as  duration  of  maintaning  the  pro¬ 
be  in  stop  position  necessary  for  conducting  reliable  measu¬ 
rements.  For  the  current  investigation  the  code  was  develo¬ 
ped  in  such  a  way  in  order  to  obtain  maximum  information 
about  flow  characteristics  in  the  vortex  region  during  mini¬ 
mum  wind  tunnel  running  time. 


1.4.  Models  of  sphere 

For  investigating  a  vortex  blunt  body  interaction  two  sphe¬ 
res  with  a  diameter  D=60  mm  were  fabricated.  One  sphere  was 
made  of  steel  and  intended  for  a  flow  visualization.  Another 
sphere  was  made  of  glass-fibre’  plastic  and  intended  for  in¬ 
vestigation  of  heat  transfer.  43  heat  flux  sensors  were  mo¬ 
unted  on  it  (Fig. 11a).  A  scheme  of  sensors  location  is  shown 
in  Fig. 11b. 
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1.5.  Heat  flux  sensors.  Method  of  heat  flux  measurement 

in  this  investigation  the  surface  thermocouples  were  used  as 
heat  flux  sensors  [18].  They  are  fabricated  from  chromel-co- 
pel  alloys.  A  diameter  of  thermocouple  wires  is  0.1  mm.  In 
the  measurement  position  they  are  being  rolled-off  into  the 
ribbons  of  0.02  mm  thickness,  0 . 2  mm  width  and  4.0  mm  long. 
The  ribbons  are  bonded  flush  with  the  model  surface.  The 
thermocouples  ends  are  brought  into  the  model  through  dril¬ 
lings  . 

Estimation  of  a  heat  flux  from  a  surface  temperature  is  a 
problem  of  the  boundary  condition  conversation.  This  problem 
may  be  considered  as  a  limited  inverse  heat  conduction  prob¬ 
lem  which  occurs  in  processing  temperature  data  measured  on 
the  body  surface.  In  such  formulation  the  inverse  heat  con- 
duction  problem  is  correct. 

As  the  model  wall  thickness  6  is  sufficiently  large,  the 
thermal  conductivity  coefficient  and  testing  time  are  small 
(i.e.  the  Fourier  number  Fo«l),  the  semi-infinite  body  model 
can  be  used  for  investigating  a  process  of  heat  transfer.  In 
this  case  the  boundary  inverse  heat  conduction  problem  may 
be  formulated  as  follows:  using  the  known  solution  T(0,t)= 
=f(X)  of  the  heat  conduction  equation 

6T(z,x)  ’  S^T  (z,x) 

— - - -  =  a  - - -  ,  'C<0,  z>0 

6t 

(here  z  Is  the  coordinate  directed  inside  the  body  perpendl- 

cular  to  the  surface) 
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with  the  boundary  conditions 

-  =  0;  T(‘^,T;)=Ti 

6z 

and  the  initial  condition 

T(z,0)=Ti 

to  determine  a  heat  flux  to  the  heated  surface 

6T(0,t) 

q(T)  =  -  I  - 

6z 

The  problem  solution  takes  the  following  form  [19]: 

\/XEp  d  X  f(t)dt  |/Icp  X  f'  (t)dt 

q(t)  = - /  =  / 

I/tT  dz  0  [/T-t  I/JT  0  1/t-t 

For  smoothing  and  functional  presentation  of  heating  curves 
K-spline  approximation  by  the  least  squares  method  is  appli¬ 
ed  [20]. 


The  thermal  activity  coefficient  value  j/Apc  incoming  into 
the  above  expression  was  specified  for  every  sensor  at  a 
special  calibration  setup  (Fig. 12a)  wherein  an  air  jet  of 
2.0  mm  diameter,  producing  a  known  heat  flux  is  directed  to 
the  sensor  perpendicularly  to  the  model  surface.  The  heat 
flux  action  period  is  approximately  the  same  as  in  the  wind 
tunnel  model  tests.  Using  the  measured  surface  temperature 
and  known  heat  flux  value  one  can  estimate  the  value  of  1/A,pc 
using  the  solution  of  one-dimensional  heat  conduction  prob- 
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lem  for  a  semi-infinite  body.  The  typical  calibration  curve 
is  shown  in  Fig. 12b. 


1.6.  Flow  visualization  methods 

For  flow  visualization  a  laser  knife  method  [15,17]  was  used 
along  with  a  shadowgraph  technique.  This  method  is  similar 
to  a  vapour  screen  method  and  allows  to  observe  a  flow  cross 
section.  For  flow  visualization  by  the  laser  knife  method 
kerosene  was  burned  up  in  the  T-33  wind  tunnel  heater.  Car¬ 
bon  dioxide  and  water  vapour  containing  in  the  combustion 
products  transit  to  a  solid  phase  during  air  expansion  in 
the  nozzle.  In  addition  there  are  particles  of  soot  in  the 
combustion  products.  Laser  with  a  power  of  6  W  and  a  wave 
length  of  5150  A  was  used  as  a  light  source.  A  plane  light 
sheet  2  mm  thickness  was  formed  by  means  of  cylindrical  op¬ 
tics.  This  light  sheet  crossed  the  gas  flow  in  the  perpendi¬ 
cular  direction  to  the  free  stream  flow.  The  flow  pattern 
pictures  were  taken  by  a  camera  located  behind  the  laser 
knife  plane.  Prior  to  the  experiment,  pictures  of  a  template 
located  in  the  laser  knife  plane  were  taken.  By  means  of 
this  template  image  scales  along  the  x  and  y-axes  were  de¬ 
termined.  During  the  experiment  a  vortex  generator  moved  re¬ 
lative  to  the  laser  knife  plane  and  in  this  way  a  flow  visu¬ 
alization  at  the  different  distances  from  the  generator  was 
carried  out. 

For  flow  visualization  on  the  sphere  surface  oil  paint  was 
iused.  The  paint  was  applied  on  the  surface  in  the  form  of 
discrete  points  [15,21].  Under  the  influence  of  a  hot  air 
flow  the  paint  is  being  softened  and  carried  away  by  the 
flow.  The  lines  appearing  on  the  model  surface  are  useful  in 
visualization  of  the  flow  direction  and  viscous  stress  va¬ 
lue  . 
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2.  Investigation  results  and  their  analysis 

2.1.  Free  stream  flow  parameters 

The  investigation  was  conducted  at  a  Mach  number  3.  The  to¬ 
tal  pressure  was  the  same  pg  =  5*10^  Pa  in  all  experiments, 
and  the  temperature  was  different:  in  one  part  of  the  expe¬ 
riment  air  was  heated  up  to  Tg  =  410  K  (by  flow  visualizati¬ 
on,  by  investigation  of  heat  transfer),  in  another  part  of 
the  experiment  the  stagnation  temperature  was  close  to  nor¬ 
mal  Tg  =  280  K  (by  measuring  pressure  and  stagnation  tempe¬ 
rature  distribution) .  In  the  first  case  the  unit  Reynolds 
number  was  Re^  =  1.1*10^  1/m,  in  the  second  case  -  Re^ 

=  3.7*10’^  1/m. 


2.2.  Flow  pattern  past  vortex  generator 
2.2.1.  Flow  pattern 

At  first  let  us  consider  a  flow  past  wing-type  generators 
and  then  past  a  vane-type  generator.  Figures  13  and  14  show 
pictures  of  vortex  cross  sections  past  generator  1  at  angles 
of  attack  CL  =  10°  and  20°.  The  image  scale  was  2:1  relative 
to  the  real  one.  Origin  of  the  x-coordinate  coincides  with 
the  rear  generator  edge.  When  analysing  the  pictures  obtai¬ 
ned  by  a  laser  knife  method  it  should  be  taken  into  conside¬ 
ration  that  in  the  flow  regions  corresponding  to  slightly 
curved  streamlines  brightness  of  the  images  is  proportional 
to  gas  density.  It  is  confirmed  by  comparison  measuring  re¬ 
sults  of  the  image  density  and  stagnation  pressure  in  the 
Prandtl-Mayer  flow  region  [17].  Due  to  this  one  can  see 
shock  waves,  rarefaction  areas,  and  a  wing  wake  in  these 
pictures.  The  results  of  digital  processing  of  the  pictures 
are  presented  in  Fig. 15.  This  processing  was  carried  out  to 
obtain  more  comprehensive  information  about  the  gas  density 
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.  fributlon  in  the  flow  field.  At  the  bottom  of  Fig. 15  a 
scale  Is  shown.  Every  gradation  on  the  scale  corres- 
to  a  certain  gas  density.  In  the  region  with  a  large 
vature  of  the  streamlines  and  particularly  in  the  vortex, 
ihese  pictures  are  not  useful  for  visualization  of  the  den- 
■tv  distribution:  due  to  the  large  curvature  of  the  stream- 
tinL  solid  and  liquid  particles  flow  away  from  the  vortex 
'ion  under  the  Influence  of  centrifugal  forces.  So  these 
Lrtices  are  seen  as  black  spots  in  the  pictures.  Similar 
paterials  about  the  vortex  formed  by  generator  2  are  presen- 

ted  in  Figs.  16-19. 

Figures  20  and  21  present  a  flow  visualization  past  genera¬ 
tor  3  which  has  the  four-vanes  form.  The  vanes  present  them¬ 
selves  rectangular  finite  span  wings.  Each  vane  produces  two 
tip  vortices.  Altogether  8  vortices  are  formed:  4  vortices 

in  the  vicinity  of  the  axis  of  symmetry  and  4  vortices  in 
the  periphery.  All  eight  vortices  and  covering  them  shock 
waves  are  clearly  seen  directly  past  the  generator  (Fig. 20, 

X  =  20  mm.  Fig. 21a).  Four  central  vortices  are  joined  into 
one  central  vortex  at  the  distance  x  -  100  mm  and  four  pe¬ 
ripheral  ones  are  retained  throughout  the  whole  region  avai- 
"'^  lable  for  observation  (x^ax  =  •  further,  for  the  in 

^vestigations  only  the  wing-type  generators  were  used  because 
they  produce  only  one  concentrated  vortex. 

Let  us  address  to  the  analysis  of  the  flow  past  the  wing-ty- 
pe  vortex  generators.  In  papers  [14,17]  a  scheme  of  the  flow 
in  the  vicinity  of  the  rectangular  wing  is  presented.  It  is 
reproduced  in  slightly  distinguished  form  in  Fig. 22a.  Ele¬ 
ments  of  this  scheme  are  seen  also  in  Fig. 14,  x=2  mm  ( 1  - 

bow  shock  wave,  2  -  internal  shock  wave,  3  -  vortex,  4 
Prandtl-Mayer  flow) .  In  the  vicinity  of  the  wing  the  vortex 
bas  a  cone  shape.  A  cone  angle  of  the  vortex  is  close  by  its 
order  of  magnitude  to  the  wing  angle  of  attack.  There  is  an 
expansion  of  the  vortex  in  the  wing  region  due  to  the  inflow 
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of  the  inviscid  gas  into  the  rotation  region.  The  flow  past 
the  wing  is  shown  in  a  schematic  presentation  in  Fig. 22b  and 
(3  past  the  wing  the  flow  expands  only  due  to  the  viscous 
dissipation.  A  vortex  expansion  angle  is  not  high  (see  be¬ 
low)  .  A  free  vortex  5  induces  a  weak  conical  compression  wa¬ 
ve  6  (Fig.  ,22c).  It  is  seen  in  the  picture  at  a  large  dis¬ 
tance  from  the  generator  (Fig. 13b,  x=135  mm;  Fig. 14b,  x-100 
mm)  Analysis  of  the  pictures  show  that  far  from  the  tip 
section  past  the  wing  a  viscous  weak  expanding  wake  7  is 
formed.  It  also  induces  a  weak  compression  wave  8  (Fig. 22b). 
Analysis  of  the  pictures  taken  at  different  distances  from 
the  generator  shows  that  the  M-lines  (Fig. 22b)  and  the 
N-lines  (Fig. 22c)  are  approximately  linear  and  located  at 
the  same  angle  relative  to  a  free  stream  flow  direction.  It 
is  close  to  the  value  of  arcsin  1/M  at  M=3.  This  confirms 
that  lines  6  and  8  observed  in  the  laser  pictures  are  weak 
compression  waves. 


2.2.2.  Vortex  dimensions 

Let  us  continue  the  analysis  of  the  laser  pictures.  The  vor¬ 
tex  axis  past  the  wing  is  curved,  the  vortex  centre  is  disp¬ 
laced  towards  the  windward  generator  surface  (in  the  positi¬ 
ve  direction  along  z-axis)  by  the  distance  up  to  0.4  b  at 
the  angle  of  attack  up  to  20-30  deg.  The  curved  portion 
length  of  the  vortex  axis  is  about  1.5b.  At  a  larger  dis¬ 
tance  from  the  generator  the  vortex  axis  is  approximately 
parallel  to  a  free  stream  flow  direction. 

Figure  23  shows  variation  of  the  vortex  diameter  versus  the 
distance  from  the  generator.  Linear  dimensions  are  related 
to  the  wing  chord.  The  vortex  has  a  spiral  shape.  It  is  cle- 
^r-ly  seen  in  Fig.  18  as  well  as  in  some  pictures  of  Fig.  24. 
The  vortex  core  has  an  oval  shape  and  its  axes  are  rotated 
t'ith  increase  of  the  distance  from  the  generator  (Fig. 24). 
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The  core  boundaries-  become  gradually,  blurred  and  the  core 
itself  becomes  not  so  black.  This  is  due  to  deceleration  of 
rotation  and  increase  of  gas  density  in  the  core.  Because  of 
insufficient  sharpness  of  the  core  boundaries  its  dimensions 
can  be  determined  only  approximately  by  means  of  the  pictu¬ 
res.  In  Fig. 23  dimensions  of  the  vortex  cross  section  in  the 
y  and  z  directions  (dy ,  )  are  shown.  At  a  =10°  and  30° 
vortex  dimensions  change  according  to  the  linear  law  throug¬ 
hout  the  range  and  at  CL  =20°  only  when  x  >  1.5  b.  The  rela¬ 
tive  vortex  diameter  value  d/b  practically  does  not  depend 
on  the  wing  dimensions:  the  experimental  points  obtained 
with  generators  1  and  2  are  on  the  general  curve  when  the 
angle  of  attack  growth.  The  vortex  diameter  d=l/2(dy+d2) 
sharply  increases  (Fig. 25). 

When  analysing  vortex  dissipation  it  is  reasonable  to  neg¬ 
lect  the  generator  dimensions.  So  in  Fig. 26  linear  dimensi¬ 
ons  are  related  to  some  conditional  reference  vortex  diame¬ 
ter  dp.  For  each  experiment  value  dp =1/2 ( dy +d2 )  was  determi¬ 
ned  at  the  distance  x=1.5  b,  where  a  transitional  curviline¬ 
ar  vortex  portion  ends  and  its  straight  portion  starts.  From 
Fig.  26,  it  can  be  seen  that  at  CL  i  20°  all  experimental  data 
regardless  to  a  generator  dimension  and  angle  of  attack  are 
described  by  a  general  function:  the  vortex  linearly  ex¬ 
pands.  An  expansion  vortex  angle  is  about  1  deg.  At  Ct=30° 
this  angle  increases  up  to  3° . 


2.2.3.  Pressure  and  stagnation  temperature  profiles 

Pressure  and  temperature  were  measured  past  generator  1 
(b=58  mm)  at  angle  of  attack  CL  =  20°  in  two  transverse  sec¬ 
tions  Xj  =2  mm  and  =110  mm  as  shown  in  Fig. 27.  The  laser 
pictures  evidence  that  the  vortex  centre  in  section  1  is 
displaced  relative  to  the  plane  z=0  towards  the  leeward  sur¬ 
face  by  13  mm.  Besides,  the  generator  is  bent  under  the  inf- 
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luence  of  aerodynamic  force  and  the  vortex  centre  is  additi¬ 
onally  displaced  in  the  direction  of  the  negative  z-axis 
approximately  by  5  mm.  Taking  into  account  the  data  obtained 
in  the  preliminary  experiments  the  probes  axes  were  set  at  a 
certain  angle  to  the  direction  of  the  y-axis  (Fig. 27a).  In 
section  one  this  angle  was  about  5°  and  in  section  two  4° . 

The  total  pressure  probe  was  rotated  in  the  plane  y=const  in 
the  angle  range  of  (p  =  ±  45°  in  increments  of  15°  .  Rotating 

was  carried  out  relative  to  the  probe  entry,  i.e.,  the  probe 
entry  position  did  not  change  during  probe  rotatings 
(Fig.27c).  At  each  position  the  probe  was  kept  for  1-2  s, 
that  provided  a  sufficient  measurement  accuracy.  The  results 
of  all  pressure  measurements  are  presented  in  Tables  1  and 
2 . 

In  Fig.  28,  as  an  example  measurements  results  of  the  stag¬ 
nation  pressure  in  section  1  and  2  at  different  distances 
from  the  vortex  centre  are  shown.  In  a  free  stream  flow 
(y'=20  mm  for  section  1  and  y'=15  mm  for  section  2)  as  ex¬ 
pected  the  probe  indications  are  approximately  symmetric  to 
the  flow  direction.  The  probe  is  practically  insensitive  to 
a  flow  deflection  up  to  the  angle  cp  ~  15°  ,  this  agrees  with 
literature  data  [22].  In  the  region  of  the  vortex  influence 
the  flow  is  swirled  in  the  clockwise  direction:  above  the 
vortex  axis  (y'=107  mm  in  section  1  and  y'=118  mm  in  section 
2)  the  flow  is  displaced  in  the  direction  of  the  positive 
2-axis  and  below  the  vortex  axis  (y'=96  mm  in  section  1  and 
y'=  93  mm  in  section  2)  in  the  opposite  direction.  Judging 
by  the  stagnation  pressure  measurements  angle  of  flow  swirl 
in  section  2  is  significantly  less  than  in  section  1. 

The  influence  of  the  rotation  angle  of  the  stagnation  tempe- 
J^ature  probe  on  its  indications  can  be  seen  in  Fig.  29  and 
tables  3  and  4.  The  stagnation  temperature  probe  indications 
practically  insensitive  to  a  flow  deflection  also  in  the 
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range  -  15°  (Fig. 29,  y'=15  mm).  The  stagnation  jbemperature 
measurements  confirm  the  presence  of  the  flow  swirl  in  a 
clockwise  direction  if  to  consider  it  in  the  flow  direction. 

each  probe  position  y=const  from  all  results  of  measuring 
pressure  or  temperature  at  different  rotation  angles  of  the 
probe  (p  the  maximum  value  was  chosen  (it  is  marked  by  an  as¬ 
terisk  in  tables  1-4).  It  was  assumed  that  this  indication 
gives  the  stagnation  pressure  value  (Ps=p, )  and  the  stagna¬ 
tion  temperature  value  (Tg=T, )  past  a  shock  wave. 

In  Fig. 30  a,  b  and  c  the  stagnation  pressure  measurement  re¬ 
sults  are  presented.  At  the  bottom  of  Fig. 30  a  and  b  ,  pic¬ 
tures  of  the  vortex  cross  section,  with  denoted  the  shock 
wave  ,  S2  and  the  vortex  boundaries  Vj  and  Vg  ,  are  pre¬ 
sented.  It  is  seen  that  the  stagnation  pressure  measurement 
results  are  adequate  to  the  flow  pattern  presented  in  the 
corresponding  picture  of  the  flow  cross  section.  In  fact,  in 
a  free  stream  flow  the  stagnation  pressure  is  approximately 
constant;  it  sharply  increases  in  the  point  S  when  crossing 
a  bow  shock  wave  of  the  generator.  Then  it  gradually  decrea¬ 
ses  while  approaching  to  the  generator,  that  agrees  with  the 
results  of  the  digital  pictures  processing.  Fig. 15.  This 
stagnation  pressure  decrease  continuously  transit  to  a  pres¬ 
sure  drop  induced  by  a  flow  swirl.  So  the  measurement  re¬ 
sults  presented  in  Fig. 30  do  not  allow  to  determine  the  vor¬ 
tex  boundary.  It  can  be  determined  by  the  value  of  a  flow 
deflection:  outside  the  vortex  dependence  of  the  pressure  p 
on  the  angle  (p  is  approximately -symmetric  to  the  free  stream 
direction.  In  section  1  the  vortex  boundaries  are  in  the  po¬ 
ints  y=95-96  mm  and  y  =  106-107  mm.  This  agrees  with  the  bo¬ 
undaries  determined  by  the  pictures  and  proves  the  use  of 
laser  pictures  for  investigating  a  vortex  development.  In 
l'lg-30c  it  is  also  can  be  seen  that  directly  past  the  gene- 
fator  (in  section  1)  the  stagnation  pressure  at  the  vortex 
axis  is  16  times  less  of  that  in  the  free  stream  flow.  In 
^oction  2  these  values  differ  only  by  2-3  times  (in  Fig. 30c 
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the  points  corresponding  to  the  nimimum  pressure  values  are 
superposed).  Hence,  with  growth  of  the  distance  from  the  ge¬ 
nerator  the  vortex  diameter  increases  slowly  but  the  flow 
swirl  decreases  sharply  and  the  vortex  core  is  filled  up 
with  a  significantly  denser  gas. 

In  Fig. 31  a  and  b  the  stagnation  temperature  measurement  re¬ 
sults  are  presented.  The  shock  wave,  as  expected,  unlike 
stagnation  pressure,  does  not  cause  increase  of  stagnation 
temperature.  However,  in  the  region  of  rotational  motion  the 
stagnation  temperature  significantly  decreases.  At  the  vor¬ 
tex  axis  this  decrease  was  18°  in  section  1  and  approximate¬ 
ly  3  in  section  2  as  compared  to  the  stagnation  temperature 
in  front  of  the  shock  wave.  The  noted  decrease  of  temperatu¬ 
res  difference  is  caused  by  deceleration  of  gas  rotation  in 
the  vortex. 

In  the  book  by  Gupta  [13]  the  measurement  results  of  air 
temperature  in  the  Ranque  tube,  described  in  paper  [23]  are 
presented.  In  this  experimental  setup  air  swirl  is  produced 
by  a  tangential  air  jet  injection  into  a  cylindrical  chamber 
at  one  its  end.  The  minimum  air  temperature  was  in  the  vici¬ 
nity  of  the  chamber  axis  at  the  same  end,  where  air  jet  in¬ 
jection  is  performed.  The  minimum  air  temperature  in  paper 
[23]  was  by  25°  lower  than  the  injection  air  temperature, 
this  agrees  with  the  results  of  the  current  experiment.  in 
the  book  [23]  gas  cooling  at  the  vortex  axis  is  explained  by 
energy  transfer  due  to  turbulent  fluctuations.  In  this  in¬ 
vestigation  the  air  flow  was  laminar.  Apparently,  air  coo¬ 
king  in  the  vortex  can  be  explained  in  the  following  way:  at 
tifst  in  the  vicinity  of  the  wing  a  vortex  with  high  circum- 
I'ential  speeds  is  formed  and  herein  speed  increases  while 
approaching  to  the  axis,  except  low  vicinity  of  the  axis, 
internal  fast  rotational  layers  due  to  gas  viscosity 
lecular  or  turbulent)  transfer  a  moment  of-  momentum  to 
external  layers  and  increase  their  circumferential  spe- 
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.  This  leads  to  decrease  of  the  internal  layers  temperatu^ 
d  increase  of  the  external  gas  layers  temperature.  As 
11  In  e  n  avers  mass  is  less  than  the  external  layers 
'  Involved  into  the  accelerated  rotation,  temperature 
“lease  of  the  internal  layers  is  essentially  higher  than 
tZ  temperature  increase  of  the  external  layers.  Local  tem- 
rature  increase  of  the  external  layer  is  seen  in  Fig^ll 
"  le  Vicinity  of  the  vortex  boundaries  V.  and  V,  and  in 
lg.31b  in  the  vicinity  of  the  vortex  boundary  V,  . 

xn  Fig  32  the  measured  stagnation  temperature  is  related 

A  the  settling  chamber  temperature  Tq  (  o 

Tortex  centres  in  sections  1  and  2  (Fig. 32,  are  conditiona  - 
ly  superposed.  The  maximum  stagnation  temperature  drop 
Ictly  past  the  generator  constitutes  only  11%.  As  a  first 
approximation  one  can  neglect  the  Influence  of  stagnation 
temperature  variation  on  gas  density  and  heat  transfer. 


2.3-  Vortex-sphere  interaction 

2.3.1.  Gas  flow  pattern  in  front  of  sphere 

^For  investigating  vortex-sphere  interaction  a  sphere  with 
the  diameter  0-60  mm  was  located  past  generator  2  (b=20  mm) 
at  the  different  distances  x,  from  its  trailing  edge.  The 
sphere  centre  was  situated  approximately  along  the  vortex 
axis.  The  investigation  was  carried  out  at  the  total  pressu¬ 
re  Po  =  5*10^  Pa  and  stagnation  temperature  Tq  -  410  K.  For 
flow  visualization  a  kerosene  heater  was  used. 

T^iftiire  of  the  sphere  in  the  absence  of  a 
Figure  33a  shows  a  picture  oi  uii  p 

vortex  generator.  This  picture  as  well  as  the  other  is  obta¬ 
ined  by  a  straight  shadowgraph  method.  Dimensions  in  e 
picture  coincide  with  the  real  ones  if  not  to  take  into 
count  small  distortions  due  to  the  light  rays  deviation. 
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p-ic  33b  a  picture  of  the  flow  past  the  generator  installed 
t  the  angle  of  attack  O  -  20»  is  presented.  Throughout  the 
“  ,ion  available  for  observation  the  flow  remains  laminar 
and  the  vortex  diameter  weakly  changes,  as  noted  in  the  pre- 

yious  chs^ptBir  • 

Figures  34  and  35  present  shadow  pictures  of  the  vortex  cord 
interaction  with  the  sphere  at  the  generator  angle  of  attack 
Ot  =  10°  and  20°  (X3=X3/D).  Analysis  of  these  pictures  shows 

the  following.  At  some  distance  in  front  of  the  sphere  a 

vortex  breakdown  occurs,  i.e.,  the  vortex  sharply  expands 

and  in  front  of  the  sphere  "a  liquid  cone"  is  formed,  which 
generates  a  cone  shock  wave.  Similar  effect  was  observed 
earlier  by  '  interaction  of  the  vortex  with  the  inlet  shock 
wave  [2],  with  the  shock  wave  in  front  of  the  cone  and  the 
barrel  shock  wave  in  the  under-expanded  jet  [12],  with  the 
shock  wave  in  front  of  two-dimensional  wedge  [24]. 


Figure  36  shows  a  supposed  flow  pattern.  The  gas  involved 
into  rotational  motion  in  the  vortex  has  a  lower  stagnation 
pressure  P3  than  the  gas  in  the  external  flow  outside  the 
vortex.  Still  a  greater  difference  is  observed  in  the  cor¬ 
responding  values  of  the  axial  momentum  component.  The  gas 
with  low  value  of  the  axial  momentum  component  can  not  pe¬ 
netrate  into  the  high  pressure  region  past  the  bow  shock  wa¬ 
ve.  so  during  flow  formation  gas  accumulation  occurs  in 
front  of  the  sphere,  followed  by  a  vortex  expansion  i.e.,  a 
vortex  breakdown.  The  flow  region  which  can  be  conditionally 
called  as  a  dead  zone  because  in  the  significant  part  of 
this  area  (in  the  vicinity  of  the  symmetry  axis)  gas  speeds 
are  low,  is  formed  in  front  of  the  sphere.  The  dead  zone  is 
a  region  of  approximately  constant  pressure.  So  it  should 
have  a  conical  shape.  In  the  point  of  intersection  of  the 
shear  layer  separating  the  dead  zone  from  the  flow  past  the 
shock  wave  with  the  sphere  surface  a  closing  shock  wave  is 
formed.  The  dead  zone  length  should  be  sufficiently  large  in 
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^(3er  the  axial  momentum  component  in  the  vortex  flow  can 
due  to  the  momentum  transfer  from  the  external  flow  to 
the  level  which  is  sufficient  for  gas  penetration  into  the 
high  pressure  region  past  the  closing  shock  wave. 

The  described  concept  is  confirmed  not  only  by  analysis  of 
the  pictures  and  also  by  results  of  calculations  based  on 
the  Euler  equations  [1].  Certainly,  such  calculations  do  not 
allow  the  liquid  cone  formation  to  be  described.  However, 
they  allow  to  determine  the  critical  relation  of  the  flow 
swirl  and  pressure  ratio  in  the  shock  wave  under  which  a 
vortex  breakdown  occurs. 

The  sphere  shock  wave-vortex  interaction  is  similar  to  the 
interaction  of  the  sphere  shock  wave  with  the  boundary  layer 
of  the  pin  installed  on  the  sphere  [25].  In  the  case  with 
the  pin  the  boundary  layer  is  separated  from  the  pin  surface 
in  front  of  the  bow  shock  wave  and  displaced  from  the  axial 
direction  to  periphery  in  front  of  the  bow  sphere  shock  wa¬ 
ve.  In  the  case  with  the  vortex  a  solid  surface  is  absent 
but  streamlines  are  also  displaced  to  periphery  in  front  of 
the  bow  shock.  Having  reached  the  sphere  surface  a  portion 
of  gas  is  deviated  inside  the  dead  zone  and  forms  a  slow  re¬ 
verse  flow  and  the  rest  portion  continues  moving  in  the  main 
direction.  On  the  sphere  surface  the  attachment  line  R  is 
formed  (Fig. 37).  In  Fig. 37  departure  from  the  axial  symmetry 
attracts  attention:  the  separation  point  S  is  displaced  re¬ 
lative  to  the  sphere  centre  0  to  the  right  side.  This  is  due 
to  eccentricity  of  the  sphere  relative  to  the  vortex  axis. 
At  the  same  time  the  attachment  line  R  has  approximately  a 
shape  of  circle  with  the  centre  at  the  sphere  axis  (in  the 
point  o) .  It  also  attracts  attention  such  a  fact  that  the 
flow  rotation  in  the  periphery  and  in  the  centre  of  the 
sphere  is  weakly  expressed. 

Analysis  of  Figs. 34  and  35  indicates  to  a  turbulent  flow  be- 
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saviour  in  the  larger  papt  of  the  dead  zone:  in  the  pictures 
(exposure  about  5  mcs)  large  vortices  are  seen  as  well  as 
compression  waves  exited  by  them.  The  separation  shock  wave 
is  an  envelope  of  these  waves.  Location  of  the  laminar-tur¬ 
bulent  transition  along  the  x-axis  changes  depending  on  the 
angle  of  attack  and  accidental  factors.  The  transition  oc¬ 
curs  either  inside  the  dead  zone  or  at  the  beginning  of  the 
dead  zone.  In  the  first  case  the  dead  zone  boundary  is  a 
broken  line  (Fig. 36  b) ,  through  to  this  point  the  boundary 
is  inclined  at  the  angle  about  6°  relative  to  the  flow  axis 
and  past  it  at  the  angle  about  20° .  In  the  second  case  the 
boundary  is  at  the  angle  of  17° .  At  a  =  20°  a  laminar-turbu¬ 
lent  transition  probably  occurs  in  the  vortex  before  the 
breakdown  point.  It  is  displayed  in  vortex  expansion.  From 
the  transition  point  a  weak  compression  wave  spreads  (see 
e  g.  Fig. 35  b) .  Evidently,  disturbances  from  the  dead  zone 
are  transferred  upstream  over  the  vortex  core.  Due  to  a 
complex  character  of  the  vortex-shock  wave  interaction  the 
interpretation  of  the  pictures  is  rather  difficult  and  the 
above  considerations  about  location  of  the  laminar  turbulent 
transition  should  be  verified. 

The  separation  shock  angle  induced  by  the  dead  zone  is  app¬ 
roximately  25°  .  It  is  close  to  the  calculated  cone  shock 
angle  [3  =  23°,  formed  in  flow  over  a  circular  cone  with  the 
apex  angle  0  =  12°  at  Mach  number  of  3. 

At  a  small  distance  between  the  sphere  and  generator  the  de¬ 
ad  zone  can  appear  between  the  leading  and  trailing  edges  of 
the  generator  as  shown  in  Fig. 34  a,  in  the  other  cases  it 
appears  past  the  vortex  generator .  In  some  cases  it  is  dif¬ 
ficult  to  determine  unambiguously  the  dead  zone  length  beca¬ 
use  of  its  complex  shape.  It  can  be  assumed  that  at  Xg/D  ) 
i-15  the  dead  zone  length  was  approximately  the  same  and 
constituted  1.05  D  in  all  experiments. 
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2.3.2.  Heat  transfer  on  the  sphere  surface 

Heat  transfer  measurements  on  the  sphere  surface  were  carri- 
ed  out  at  X3  .  70  mm  (X3/D=1.17)  and  three  angles  of  attack 
of  the  vortex  generator  2:  a  =10» ,  20»  and  30« .  The  measure¬ 
ments  on  the  sphere  surface  were  also  carried  out  without  a 
vortex  generator.  The  experements  were  performed  at  Tj,  =  400 
K.  An  electrical  air  heater  was  used. 

The  sphere  was  located  in  the  wind  tunnel  test  section  in 
such  a  way  that  the  generatrix  with  thermocouples  1-34  was 
in  the  plane  z=0,  and  the  generatrix  with  thermocouples 
35-43  -  in  the  plane  y=0.  The  experiments  were  performed  in 
the  following  way.  During  the  period  of  setting  a  flow  regi¬ 
me  in  the  wind  tunnel  the  model  was  in  the  thermostatic 
chamber  and  was  separated  from  the  flow  by  a  special  flap. 
Airflow  over  the  model  surface  was  performed  by  a  cool  air 
through  the  pipe,  situated  near  the  frontal  surface.  After 
setting  a  prescripted  flow  regime  the  model  was  injected  in¬ 
to  the  flow  by  means  of  the  fast-acting  mechanism.  The  model 
has  been  in  the  flow  for  3  seconds.  During  this  period  re- 
'  cording  of  thermocouples  signals  was  carried  out  with  samp¬ 
ling  rate  of  100  1/s,  then  the  model  was  took  out  from  the 
flow  and  the  wind  tunnel  operation  was  stopped. 


The  results  of  the  measurements  are  presented  in  Tables  5  8. 
The  following  parameters  are  presented  in  the  tables:  the 
generator  angle  of  attack  tt,  the  Mach  number  M,  the  Reynolds 
number  Re,  the  stagnation  temperature  Tq ,  the  total  pressure 
Po  •  The  time  T  =  1  s  is  also  presented  there,  for  this  time 
the  obtained  data  processing  was  being  carried  out.  The  cho¬ 
sen  time  is  essentially  larger  than  the  period  of  injection 
of  the  model  into  the  flow  and  it  is  quite  small  to  avoid  of 
noticeable  heat  spreading  in  the  region  of  the  maximum  heat 
transfer  on  the  sphere  surface.  In  the  first  column  of  the 
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tables  the  thermocouple  number  is  pointed  out,  in  the  second 
one  -  its  coordinate  s/R,  further  -  the  actual  surface  tem¬ 
perature  T  at  X  =  1  s,  the  smoothed  temperature  Tgmo ^  the 
jieat  flux  q,  the  calculated  value  of  a  heat  flux  in  the 
stagnation  sphere  point  qo  at  the  temperature  T  and  the  re¬ 
lation  q/qo  • 

The  heat  flux  in  the  stagnation  sphere  point  was  calculated 
from  the  Fay-Riddell  formula  [26]  which  for  the  perfect  gas 
takes  form 

qg  =  0.94(p,,/i^,)°-^  (pg/lg  (ig -cT^,  )R"° (6u/6s)° 

In  hypersonic  flow  over  the  sphere  for  calculation  of  the 
velocity  gradient  in  the  stagnation  point  we  can  use  the 
formula  obtained  under  condition  of  the  Newton  pressure 
distribution  on  the  surface 

6u/5s  =  l/R(2(Ps-p^)/Ps (1) 

In  the  current  paper  for  improving  accuracy  of  the  pressure 
gradient  calculation  at  M=3,  calculation  results  for  invis- 
cous  flow  over  the  sphere  were  used  [27].  The  difference  in 
the  value  5u/5s  as  compared  to  formula  (1)  was  7.5  %. 

Distribution  of  the  relative  heat  flux  q/qg  over  the  sphere 
surface  in  the  case  with  a  laminar  flow  was  calculated  by 
the  equation  [28] 

q/qo  =  qo  +  giCos((i)n)  +  g2Cos(2(DTi)  +  gs  Cos  (30)11) 

Here  Ti  =8e’^(s/s,  ),  s,  is  a  sonic  point  coordinate;  n=0.5; 
W=i;  gg  =0.452;  g^  =0.499;  gg  =0.049;  g3=0;  36=1.4.  Due  to  acco¬ 
unt  for  the  influence  of  the  Mach  number  on  the  sonic  point 
location  this  formula  gives  a  higher  accuracy  than  a  simple 
^nd  often  used  equation 
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q/qg  =  0.55  +  0 . 45Cos ( 2s/R ) 

In  this  paper  a  sonic  point  coordinate  was  determined  by  the 
tables  [27]. 

In  Fig. 38  calculated  and  experimental  heat  flux  distributi¬ 
ons  over  the  sphere  surface  in  the  absence  of  a  vortex  gene¬ 
rator  are  compared.  The  heat  flux  values  are  presented  for 
the  surface  temperature  T,^=290K.  It  is  seen  that  under  these 
conditions  a  laminar  flow  over  the  whole  sphere  surfaces  was 
retained.  Some  divergence  of  the  heat  flux  values  at  the 
differente  generatrics  of  the  sphere  as  well  as  the  diffe¬ 
rence  of  a  calculated  and  experimental  values  in  the  vicini¬ 
ty  of  the  stagnation  point  are  due  to  the  experimental  er¬ 
ror.  Furthermore,  excess  of  the  experimental  value  over  the 
calculated  one  in  the  vicinity  of  the  stagnation  point  can 
be  associated  with  incoming  flow  turbulence:  according  to 
the  measurements  carried  out  by  a  laser  Doppler  Velocity  Me¬ 
asurement  System  (DLA)  the  flow  turbulence  in  the  T-33  wind 
tunnel  constitutes  approximately  1%. 

Figure  39  shows  the  influence  of  the  vortex  cord-bow 
wave  interaction  on  the  heat  flux  distribution  over  the 
sphere  surface.  In  Fig. 39  both  calculated  and  experimental 
heat  flux  values  are  related  to  the  calculated  heat  flux  in 
the  stagnation  point  of  the  sphere.  The  vortex  causes  a  qua¬ 
litative  variation  of  the  heat  f.lux  distribution:  in  the  vi¬ 
cinity  of  the  sphere  axis  the  heat  flux  is  substantially 
less  than  on  the  isolated  sphere.  In  the  vicinity  of  the  at¬ 
tachment  line  R  the  heat  flux  in  the  presence  of  a  vortex 
essentially  increases  not  only  as  compared  to  the  heat  flux 
It  the  same  point  and  in  comparison  with  the  heat  flux  in 
the  stagnation  point  of  the  isolated  sphere  (Fig.  39b,  CL 
•  Such  a  significant  heat  flux  increase  can  be  explai¬ 
ned  by  two  things: 
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Dsubsequent  gas  compression  in  the  oblique  and  normal  shock 
waves  leads  to  a  greater  pressure  increase  at  the  attachment 
line  R  than  the  gas  compression  in  the  normal  shock  in  the 
vicinity  of  the  stagnation  point  of  the  isolated  sphere; 
2)vortex-bow  wave  interaction  leads  to  a  flow  turbulation  in 
the  vicinity  of  the  sphere.  The  flow  turbulence  also  leads 
to  a  significant  increase  of  the  heat  transfer  downstream  of 
the  line  R  on  the  cylindrical  part  of  the  model  (compare 
q/qo  at  s/R~2.5,  Fig. 38b  and  39b).  As  a  result  of  the  vortex 
influence,  an  integral  heat  flux  to  the  blunt  body  surface 
increases  as  well. 

The  heat  flux  distribution  on  the  sphere  interacting  with  a 
vortex  is  similar  to  the  heat  flux  distribution  on  the  spi¬ 
ked  sphere  [25]:  in  case  two  the  heat  flux  in  the  vicinity 
of  the  symmetry  line  also  decreases. and  increases  at  the  at¬ 
tachment  line  and  behind  it. 

From  Fig. 39  it  is  seen  that  the  heat  flux  distribution  over 
the  sphere  surface  interacting  with  a  vortex  was  asymmetri¬ 
cal  at  all  investigated  values  of  the  generator  angles  of 
attack.  Probably,  it  is  associated  with  that  fact  that  the 
vortex  and  sphere  are  not  coaxial.  Judging  by  Fig. 39  at  all 
angles  of  attack  the  vortex  was  essentially  displaced  rela¬ 
tive  to  the  sphere  axis  in  the  vertical  plane:  at  ”K  =  180°  the 
heat  flux  maximum  is  substantially  less  and  displaced  by  a 
greater  distance  from  the  axis  of  symmetry  than  at  ')(=0  (Fig. 
39).  It  can  be  the  consequence  pf  that  fact  that  the  vortex 
axis  goes  below  the  sphere  axis. 

In  Fig. 40  the  heat  flux  distributions  at  different  angles  of 
attack  of  the  generator  are  compared  for  each  sphere  gene¬ 
ratrix.  The  maximum  relative  heat  flux  values  qm/qo  versus 
the  generator  angle  of  attack  are  presented  in  Fig.  41.  From 
Figs.  40  and  41  it  is  seen  that  qm/qo  changes  against  the 
angle  in  nonmonotonous  manner:  at  (X  <  20°  the  value  qm/qo 
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,o«s  with  Increase  ct  the  angle  and  at  a  higher  angle  of 
\tacK  it  decreases.  However,  this  effect  can  not  be  unambi- 
Lously  associated  either  with  the  flow  swirl  increase  or 
the  vortex  diameter  increase  against  the  angle  of  attack. 
Apparently,  the  value  q„/qo  was  greatly  influenced  by  the 
rtex  axis  displacement  relative  to  the  sphere  axis.  At  a 
-20”  the  eccentricity  was  the  least  one  (coincidence  of  heat 
nux  distributions  at  K-0  and  T(=90“  is  an  evidence  of  this), 
probably,  it  was  the  reason  that  the  maximum  relative  heat 
flux  value  was  obtained  at  (X=20  . 


Por  quantitative  comparison 
interacting  with  the  vortex 
on  the  ordinate  axis  of  Fig. 
paper  by  Crawford  described 
mum  values  obtained  in  this 
ford  results . 


of  heat  transfer  on  the  sphere 
and  on  the  sphere  with  the  pin 
41,  the  maximum  value  from  the 
in  [25]  is  presented.  The  maxi¬ 
investigation  are  close  to  Craw- 


Thus,  the  measurements  have  shown  that  due  to  the  vortex 
influence  the  flow  over  the  blunt  body  is  being  reconstruc¬ 
ted  and  on  its  surface  a  significant  heat  transfer  increase 
is  observed,  that  is  very  important  for  practical  applicati- 
ons . 
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Nomenclature  for  the  part  B 


a 

Cp 

Cv 

g 

h 

ho 

M 


R 

S 

T 


=  of  sound, 

=  specific  heat  at  constant  pressure, 

=  specific  heat  at  constant  volume, 

=  velocity  circulation  in  axisymmetric  flow, 

=  Yp/(  ("^-D/P)  ^  enthalpy  per  unit  mass, 

=  h+(u2+v2+w2 )/2,  total  enthalpy  per  unit  mass 

=  Mach  number, 

=  pressure, 

=  radial  coordinate, 

=  u^  '  Reynolds  number, 

=  entropy  per  unit  mass, 

=  absolute  temperature. 


=  axial,  radial  and  azimuthal  components  of 

X  =  coordinate  along  free  stream  direction, 
y  =  r^/2, 

*)(  =  Cp/c^,  ratio  of  specific  heats, 

5  =  vacuum  core  scale, 

5  =  effective  vortex  thicknesses 

1\  =  r^/(4x)  , 

jl  =  viscosity  coefficient, 

jig  =  second  viscosity  coefficient, 

^  ^  (i-2_i)/2, 

p  =  density, 

6  =  Prandtl  number , 

=  time, 

*P  =  swirl  angle, 

^  =  stream  function, 


velocity , 


g.h, 

P.T  =  circulation,  enthalpy,  pressure 
by  external  values  at  r 


,  temperature  scaled 
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„  =  total  temperature. 

lo 

subscripts 

CO  =  free  stream  conditions, 
j  =  first  term  in  asymptotic  expansion, 

2  =  second  term  in  asymptotic  expansion, 

(j  =  dimensional  values, 

e  =  external  flow, 

i  =  initial  conditions, 

s  =  stagnation  point  past  normal  shock  wave. 
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B 


.  walitical  ana  numerical  study  cf  vorten  development 


1.  problem  formulation 

„ain  theoretical  approaches  used  for  description  of  the 
"  “  development  and  hreahdown  in  incompressihle  fluid 

Ho  41)  are  generalited  here  to  include  compressibility.  The 
s  -  ylindrical  approximation  is  derived  as  asymp to 
rlacl  for  the  slender,  smoothly  developing  vortex.  Small 
turbations  of  the  cylindrical  flow  are  considered  to  ge 
Terau"  a  well-Known  classification  of  the  vortex  flows  on 
subcritical  and  supercritical  [30]. 

1.1.  Governing  equations 

•will,  of  oerfect  gas  having  axial, 

Axisymmetric  compressible  flow  o 

radial  and  azimuthal  components  of  velo  y  .^^^tions 

The  viscosity  and  the  bulk  viscosity  ^2  ^re 

Of  temperature,  the  Prandtl  number  6  and  K-C,/C.  are  o 

tants  The  Navier-Stokes  equations  for  axial,  radial 

,  and  the  mass-conservation  equation  are 

azimuthal  momentums  and  rne 

represented  correspondently  as. 


du 


dX 


6p  1 

—  +  /iV^u  +  -  (M+2M2)  "" 

6x  3  6x 


(1 .1  -1) 


5m 

6u  2 

6v  6u 

2  S|l2 

+  - 

(2 - -  D)  +  — 

(—  +  — ) 

+ 

5x 

6x  3  Sr 

6x  Sr 

3  Sx 

D  , 
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dv 

W 

6p 

V 

1 

8d 

f - - 

_  N  ^ 

-  — 

+  M 

(V^v 

-  — ) 

+  — 

(M+ 

2M2) 

p  ( 

dt 

/ 

r 

5r 

r2 

3 

8r 

8m 

6v 

8u 

8m 

8w 

2 

2  8)12 

+  —  < 

'  —  + 

— ) 

+  - 

(2— 

-  -D) 

+ 

—  — 

-  D 

8x 

5x 

8r 

6r 

8r 

3 

3  8r 

dw 

w 

w 

8)1 

8w 

8)1 

8w 

w 

p  (—  ^ 
dt 

r  V - ) 

r 

=  M 

-  — ) 
r2 

+  — 

8x 

—  + 

8x 

5r 

{— 

6r 

-  -) 

r 

+ 


(1.1.2) 


,  (1.1.3) 


ii 


(1.1-4) 

—  -1-  pD  =  0  , 

dX 

where  x,r  -  cylindrical  axial  and  radial  coordinates;  u,v,w 
-  axial,  radial  and  azimuthal  components  of  velocity; 

p-pressure ;  p  ““  density, 

d  6  6  6 

—  =  —  +  u  —  +  V  , 

dt  8x  6r 

6u  1  6 

D  =  div  V  =  +  “  ( vr )  , 

6x  r  6r 

6^  8^  1  5 

V2  =  -  +  -  + - . 

8x2  5j-2  j.  5r- 


1 

k 


il 


Ii 

tl 


I 
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we  reproduce  here  also  equations  for  enthalpy  h-'lfp/((')(  l)p) 
and  pressure; 

dh  dp  P- 

p  —  =  —  +  div(  -  Vh)  +  N  ,  (1-1-5) 

dt  dt  6 


dp 

—  =  -  ')(pD  + 
dt 


('^-l)div(  -  Vh)  +  ('K-l)N 
6 


5u  6v  6w 

N  =  JLl  [2(  — +  2(  — +  2  —  +  (  — )^  + 
5x  6r  5x 


(1.1.6) 


5v  6u  6w  w  2 

+  (—  +  — )2  +  (—  -  -)2]  -  -  (M-M2)d"  . 

5x  6r  6r  r  3 

The  equations  (l.l.l)-{l-l-4)  together  with  any  of  equations 

(1.1.5),  (1.1.6)  and  with  stated  functions  /1(h),  /12(h)  compo¬ 

se  the  complete  system  which  is  to  be  analyzed. 

1.2  Inviscid  axisymmetric  steady  flow 

The  inviscid  counterparts  of  the  equations  considered  above 
are  obtained  setting  /1=0,  /Ig  =0 .  Then,  the  well-known  consequ¬ 
ences  of  the  equations  are: 

6hQ  Shg 

u  -  +  V  -  =  0  , 

6x  5r 
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6s  6s 

u  —  +  V  =  0  , 

6x  6r 


(1.2.1) 


6g  6g 

u  —  +  V  —  =  0  , 
6x  6r 


where  ho =h+(u2 +v2 +w2 )/2  -  total  enthalpy;  S=So +0^ ln(p/p^ )  - 
entropy;  g=wr  -  circulation.  The  mass-conservation  equation 
(1.1.4)  is  integrated  by  introducing  of  the  stream  function 


'll): 


S')]) 

pur  =  —  , 
6r 


(1.2.2) 


pvr  =  -  —  . 

6x 


The  equations  (1.2.1)-(1.2.2)  demonstrate  that  the  functions 
ho,  S,  g  are  convected  along  streamsurf aces .  Therefore, 
three  integrals  exist  for  the  inviscid  axisymmetric  steady 
flow:  ho=ho('i|)),  s=s(i|)),  g=g(il)). 

We  notice,  that  entropy  is  always  convected  together  with 
gas  particles,  total  enthalpy  is  convected  along  the  stream¬ 
lines  only  in  steady  flows.  The  circulation  in  arbitrary 
three-dimensional  flow  is  described  by  following  equation 


dr 

—  =  f  (VTXVs)n  dA  , 
dX 


1 
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Where  r  .  -  circulation  along  a  closed  curve  formed  by  moving 
articles,  T-temperature ,  dA  -  area  element  of  the  surface 
based  on  this  closed  curve,  TT  -  unit  vector  normal  to  this 
element.  Therefore,  in  contrast  to  incompressible  flow,  the 
circulation  is  conserved,  and  the  conservation  equation 
(1  2.1)  is  ^  consequence  of  axisymmetry . 

using  the  momentum  equation  represented  in  the  Lamb's  form 

V2  ^-4  1 

V( — )  +  (rot  V)XV  =  -  -  Vp  , 

2  P 

1 

and  the  thermodynamic  correlation  tVs  =  Vh  -  -  Vp  ,  the  fol¬ 
lowing  relation  is  found:  P 


VXrot  V  =  Vhg  -  tVs  , 


The  right-hand  side  vector  in  (1.2.3)  is  normal  to  the  stre¬ 
am  surfaces.  Therefore,  the  vorticity  Q  =  rot  V  has  zero 
component  in  this  normal  direction,  its  component  along  the 
meridian  section  of  the  stream  surface  -  pq5  ( wr  ) /6'\|),  q 
(u^  +  and  azimuthal  component  Qy=  6v/6x  -  6u/6r.  The 

velocity  vector  is  decomposed  as  \(j-=  g>  Vjj  -0,  ■  Using 

these  coordinates,  (1.2.3)  is  represented  as. 


6v 

5u 

dS 

dho 

1  dg^ 

—  - 

—  =  pr  (T  — 

- 

+ 

6x 

6r 

dlj) 

d'lj) 

2r^  dlj) 

We  use  the  following  identities  to  transform  the  left-hand 
part  of  (1.2.4)  into  adapted  form: 
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.P 

S(ll))  =  So  +  Cvln-  , 

P"' 

P 

^2  =  “X-  ,  T  =  aV[(^-l)Cp3  • 

P 


Therefore,  the  invisoid  axisymmetrlc  compressible  steady 
flow  is  described  by  single  quaslllnear  second-order  equati¬ 
on  (1  2.5).  As  far  as  the  shocks  may  be  formed  in  the  flow, 
the  equations  must  be  supplemented  by  the  conditions  on  the 
shock  we  do  not  concern  this  question  here  and  only  note, 

that  the  functions  h„  (i|))  ,  gfip)  do  not  change  through  the 

Shock,  while  S{^)  has  a  jump. 


The  simplist  application  of  equations  (1.2.5),  (1.2.2) 
(1.2.6)  is  the  investigation  of  the  cylindrical  vortex  res 
ponce  on  the  changes  of  the  external  pressure. 


Let  the  axial  scale  is  large  compared  with  the 
and  V  «  u.  The  small  terms  are  discarded  in 
the  second-order  ordinary  differential  equation 


radial  one 
(1.2.5)  and 
is  obtained; 


d  1  dlj)  dho  i  PH*) 

r  — ( - )  =p^r2  [ - (I]))  - - 

dr  r  dr  dH* 

(1.2.7) 

dS('ll))  pug^(l|)) 

-  T  - - (i  +  ('K-i)  — )]  +  ” 

,2 


Let  the  initial  profiles  Ui(r)>0,  gi  (r) ,  Pi(r),  and  the 
pressure  at  ifinity  (cc)  are  determined  at  some  section . 

The  pressure  and  stream-function  are  found  by  integration: 


( 
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00  p^g^2 

Pi  (r)  =  Pi  (~)  -  /  - 


^  (r)  =  /  piUirdr  . 

0 

Then,  the  inverse  function  r{'^)  is  obtained  and  the  functi¬ 
ons 

gC'ip)  =  91  ' 

Pi 

,  (1.2.9) 

S(l|))  =  So  +  Cyln 

p/  r=r(il)) 

'y  Pi  ui^  91^(11)) 

ho(il))  =  ( - ^  ^ 

^-1  p,  2  2r2  r=r(ip) 

are  determined. 

The  profiles  p,  u,  p,  T  are  found  at  v=0  from  (1.2.2), 
(1.2.6),  like  it  was  mentioned  above,  and  are  the  functions 

of  ij),  r,  d'\|J/dr . 

The  profiles  g,  S,  hg  are  considered  tending  to  constants  as 
ijJ  00,  i.e.  the  vortex  is  surrounded  by  potential  flow.  Let 
the  external  pressure  changes  from  Pi (“)  to  p(“).  Then,  the 
velocity  and  density  at  infinity  are  obtained  from 
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K  Pi(co)r  P(“) 

(CO)  =  +  2  1- 

K-l  pi  (00)  L  Pi  (00) 

-  P(“) 

p(co)  =  p^  (00)  - 

Pi  (“) 

The  responce  of  the  vortex  on  this  external  pressure  change 
is  found  by  solving  of  (1.2.7)  at  boundary  conditions 

1  Sip 

l|j(0)  =0, - ^  p(oo)u(oo)  at  r  oo  . 

r  6r 

It  is  easy  to  show  that  the  profile  'i|)(r)  found  from  (1.2.8) 
is  a  solution  of  (1.2.7)  for  boundary  conditions 

1 

l|)(0)  =0, - ^  pi(oo)Ui(oo)  at  r  ^  00  . 

r  5r 

Do  other  solutions  of  (1.2.7)  exist  at  the  same  boundary 
conditions  ?  Such  solutions  were  found  in  incompressible 
flows  and  were  called  "conjugate"  to  the  basic  one  [30] .  The 
existence  of  conjugate  flows  allowes  spontaneous  jumps  bet¬ 
ween  them  which  were  treated  as  vortex  breakdown. 

1.3  Small  perturbations  of  the  cylindrical  flow 

The  analysis  of  steady  small  perturbations  of  the  cylindri¬ 
cal  vortex  performed  for  incompressible  flow  resulted  in  an 
important  division  of  the  vortex  flow  on  subcritical  and 
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supercritical.  We  generalize  this  analysis  for  the  compres 
sible  vortex. 

The  governing  equation  may  be  derived  from  (1.2.5).  However, 
it  is  useful  to  perform  the  derivation  directly  from  the  Eu¬ 
ler  equations: 

6u  6u  6p 

p  (u —  +  V — )  =  -  —  , 


6x  6r 


5x 


6g 

6g 

u —  + 

V — 

=  0  , 

5x 

6r 

6p 

6p 

6u 

U -  + 

V - 

5x 

6r 

5x 

6(pu) 

1 

6(pvr) 

■)  , 


6x 


r  5r 


(1.3.1) 


6(vr)  6(vr)  g^+v^r^  6p 

p  (u -  +  V -  -  - )  =  -r — 


6x 


6r 


6r 


The  small  perturbations  of  the  basic  cylindrical  flow  are 
considered : 

«  =  Uq  (y)  +  (x,y)  +  . . . 

g  =  go  (y)  +  ^gi  (x,y)  +  •  •  •  . 


+  ,  .  .  , 


p  =  Po  (y) 


Epi (x,y) 


(1.3.2) 
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p  =  Po (y)  +  ^Pi  (x,y)  +  .  .  . 

V  =  (X,y)  +  .  .  . 

where  V=vr ,  y=r^/2  and  the  base  cylindrical  profiles 

are  correlated  by  the  equation 

dPo  Po  9o 


dy  4y^ 

The  equations  for  perturbations  are: 

6ui  6pi 

Po  (^0 -  +  Vi  ^0  '  )  =  -  -  , 

6x  6x 


(1.3.3) 


6gi 

Uq  +  Vj  qQ  '  —  0  , 

6x 

Su^  6pi  Svj 

Po -  +  Uo -  +  Po  -  +  Vjpo'  =  0  ,  (1.3.4) 

6x  6x  6y 

6pi  6ui 

Uo -  +  VlPo'  =  -ypo( -  +  - )  , 

6x  6x  6y 

90^  gggi  6pi 


-pi -  +  Po(Uo -  -  - )  =  -2y -  , 

2y  6x  y  6y 

where  the  prime  means  the  y-derivative . 

The  axial  derivatives  of  the  perturbations  are  found  from 
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v’the  former  four  equations  of  (1.3.4): 

1  6pi  Po  ' 

-  =  ( -  f  +  f’)/(l-Mo2)  , 

PqUo^  6x  Ypo 

1  6ui  Uq  '  Po  ' 

- -  -  -  f  -  ( -  f  +  f')/(l-Mo2)  , 

Uo  6x  Uo  ^(Po 

6gi 

—  =  -  go  'f  . 

5x 


1  5pi  Po  '  Po  ’ 

- =  -  -  f-f  +  ( -  f  +  f')/(l-Mo2)  , 

Po  5x  Po  "KPo 

where  f  =  /Uq  ,  Mq  ^  =  Po^o^/CKPo)  • 

Using  (1.3.5)  and  the  last  equation  of  (1.3.4),  the 
wing  second-order  linear  equation  is  obtained: 

PoUo^f '  '  poUo^  5^f  Po  ’  '  2  Po  '2 

( - )  + - +  [( - )  +  -  Po  ' - 

2y  6x2  y  KPod-M, 


The  ordinary  differential  equation  is  found  for  the 
bations  having  a  form  f  =  exp((Xx)  F(y)  : 


(1.3.5) 


follo- 

- ]f  =0. 

(1.3.6) 

pertur- 


-F’  )  '  +  t' 


+  (■ 


l-Mo' 


KPo(1-Mo") 


(1.3.7) 


r^oint  M  =1  is  a  singular  point  of  equation 

we  must  Investigate  this  singularity,  especially 

^  4-«v.«,  /1-M  is  a  denominator 

taking  into  account  that  the  term  (1  Mo  )  is 

in  (1.3.5). 

ikt  Ho-l  at  y=y„  and  8=y-y„ •  The  equation  (1.3.7)  is  repre- 
sented  near  d=0  as 


(H.a9*...)F"  -  (l+B9t...)F'  *  (CtD8t...)F 


=  0  , 


B  =  -  A 


-  (Mo")'  ^ 


(1.3.8) 


c  =  - 


D  =  (')(-l)C^  -  C(-  A  -  (Mo^  )  '  )  7 


v/here  the  coefficients  A,B,C,D  are  ca 


Iculated  at  d=0. 


According  to  the  method  of  Frobenius  ,  two  linearly  indepen 
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^ent  solutions  of  (1.3.7)  near  ^  =  0  have  the  form: 

00 

Fi  =  ^  an  a"  , 

n  =  0 

00 

Fg  =  QFi  Inia  I  +  X  bnd"  , 

n  =  0  • 

where  coefficients  an,  bn,  Q  are  determined  from  (1.3.7): 

1 

ag  =  1,  =  -  -(2(A-B)+C),  ...  , 

3 

bg  =  1 ,  bj  =  C,  ...  , 

1 

Q  =  -  -[D+(B-C)bi  ]  . 

2 

In  general  case,  the  solution  has  infinite  second  derivative 
at  a=0.  However,  direct  calculation  showed  that  Q=0,  and  all 
derivatives  are  finite.  Moreover,  the  generally  singular 
term  (fpg  '/KPo ' )/( I'^g^ )  is  regular.  , 

An  arbitrary  jump  of  the  second  derivative  is  allowed  by  the 
equation  (1.3.7)  at  y-y^  •  But  in  this  case,  the  functions 
'Ji  ,  Pi  ,  Pi  would  have  a  jump  too,  and  we  exclude  this  possi¬ 
bility.  This  requirement  makes  the  problem  closed. 

The  another  question  is  the  asymptotic  behaviour  of  the  so¬ 
lutions  of  (1.3.7)  at  y  00.  For  this  purpose  we  must  know 
’’'ore  details  about  structure  of  the  basic  cylindrical  flow. 

consider  the  vortex  surrounded  by  potential  flow  with 
l^o=const,  S=const,  g=const.  The  deflexion  of  these  functions 
Ifom  constants  is  caused  inside  the  vortex  by  diffusion  and 
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is  usually  exponentially  small  far  from  the-  vortex  core.  The 
potential  vortex  has  its  own  dimensional  scale 


a«^  2  ^ 

where  a^  -  the  sound  velocity  far  from  the  core.  The  soluti¬ 
on  is  "easily  found  from  (1.3.3)  and  the  conditions 
=const ,  S=const,  g=const: 


Uo  = 


nj  =  nj/{i  -  y*/y)  . 


io  y* 


(1.3.9) 


Po 


y« 

1  -  - 

y 


■x/df-i) 


Po 


y. 

1  -  - 

y 


i/(')(-i) 


This  solution  represents  a  vortex  having  the  vacuum  core  at 
y=y,  ,  Mq  as  y  y*  •  The  solution  has  not  counterparts  in 
incompressible  fluid  (a  rather  remote  analogue  is  the  vortex 
with  cavitation).  The  expansion  of  the  solution( 1 . 3 . 9 )  at  y 
»  y*  is  the  general  representation  of  the  flow  around  a 
compressible  vortex.  Using  this  representation,  the  asympto¬ 
tic  behaviour  of  the  solutions  of  (1.3.7)  at  y  “  is  found; 
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F  =  (2y)^^^exp(6l/2y)  [1  -  C(2y)'^^^  +  0(y'M] 

P  =  ±  (M«.^  -  1))^''^  f  (1.3.10) 
-  3  2 

c  =  -  —  +  Py*  ( — :  -  - ) 

L  sp  y-i  1-n.J- 

Two  signs  in  (1.3.10)  correspond  to  two  different  linearly 
independent  solutions.  The  positive  eigenvalues  OC^>0  descri¬ 
be  the  upstream  influence  of  the  profile  perturbed  at  some 
position,  while  0?  <0  describe  the  standing  waves.  In  two  ca¬ 
ses:  (X^>0,  Mao>l  and  0?  <0 ,  Moo  <1 ,  one  of  two  functions  incre¬ 
ases  exponentially  as  y^a>  and  obviously  must  be  excluded. 

The  following  physical  examination  is  usefull  for  understan¬ 
ding  of  possible  solutions  of  (1.3.7).  We  consider  the  vor¬ 
tex  confined  into  finite  cylinder  (0<y<yR).  The  boundary 
conditions  for  (1.3.7)  in  this  case  are:  F(0)=0,  F(yp)=0. 
Only  monotonous  profiles  of  Mg  are  considered.  The  asserti¬ 
ons  are  not  rigorously  proved  and  are  obtained  from  asympto¬ 
tic  analysis  of  the  solutions  for  |CX^|>1. 

1.  In  the  case  Mo^<l  the  infifnite  discrete  spectrum  exists 
such  that  (Xg  ^  <0Cj  ^  <(I2  ^  <  ...  .  The  positive  OCg^>0  are  conden¬ 
sed  into  continuous  spectrum  as  y^  ^  oo.  In  this  limit  only  a 
finite  number  of  discrete  eigenvalues  io  may  exist.  If  at 
least  one  OC^  <0  exist  -  the  flow  supports  standing  waves  and 
is  called  "subcritical" ,  in  opposite  case  -  "supercritical". 
This  is  direct  generalization  of  the  classification  adopted 
for  incompressible  flows. 

2-  In  the  case  Mg^>l  the  infifnite  discrete  spectrum  is  or- 
^®5^ed  as  CCg  ^  XX^  ^  >(X2  ^  >  ...  .  The  standing  waves  are  always 
supported  in  the  confined  supersonic  flow.  They  are  caused 
simply  by  reflection  from  the  boundaries.  For  physical  rea- 
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sons  all  CC^  are  negative  (however,  we  did  not  prove  it  di- 
rectly  for  the  equation  (1.3.7)).  All  negative  eigenvalues 
Qj;e  condensed  into  continuous  spectrum  as  yjj  ^  oo.  is  there  a 
positive  eigenvalue  possible  ? 

3.  In  the  case  Moo  Mo  (0)<1,  the  discrete  spectrum  covers 
whole  real  axis:  .  .  .  <(X_  ^  <(Xo  ^  <(Xi  ^  <...  .  Its  negative  part 
is  condensed  into  continuous  spectrum  as  y^  «».  The  infini¬ 
te  discrete  positive  spectrum  always  exists  for  unbounded 
flow.  From  this  spectrum  we  can  determine  how  the  profile 
perturbations  propagate  upstream. 

4.  In  the  case  <1,  Mq  (0)>1  the  discrete  spectrum  covers 
whole  real  axis  too.  Its  positive  part  is  condensed,  however 
infinite  discrete  negative  spectrum  remains  as  yj^  co.  such 
unbounded  flow  always  supports  standing  waves  attributed  to 
the  reflection  inside  the  inner  supersonic  part. 

Therefore,  two  problem  statements  concerning  discrete  spect¬ 
ra  in  unbounded  flow  are  formulated:  if  Ma,<l  -  to  find  all 
ft^<0,  and  if  M^>1  -  to  find  all  0? >0 ,  such  that  the  equation 
(1.3.7)  and  boundary  conditions  F(0)  =  0,  F'(oo)=o  are  satis¬ 
fied. 

The  most  important  finding  possible  during  the  discrete 
spectra  analysis  is  the  detection  of  an  eigenvalue  which 
tends  to  zero  as  some  profile  parameter  changes.  The  exis¬ 
tence  of  the  eigenvalue  close  to  zero  means  that  the  pertur¬ 
bations  propagate  far  upstream. 

The  asymptotic  behaviour  of  the  solutions  at  (X^=0,  y  oo 
differs  from  (1.3.10).  Two  linearly  independent  solutions 
have  the  form: 
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ai  1 


P  =  1  +  +  ) 

y  y 

Iny 

=  y  +  bi  Iny  +  0(  ) 


ai  - 


(1.3.11) 


y*  ' 


bi  = 


M. 


•2ai  -  (■ 


1-Mc 


■)y. 


^-1 


1.4  Quasi-cylindrical  approximation 

he  quasi-cyllndrlcal  approximation  is  an  asymptotic  appro- 
ch  for  description  of  the  slender  vortex  smoothly  regrou 
dng  under  the  action  of  viscosity  or  pressure  gradient, 
s  similar  in  its  nature  to  the  boundary-layer  approach. 

.et  the  vortex  has  a  radial  scale  5  in  its  initial  section, 
londimensional  variables  are  introduced  by  the  scaling: 
ra=6r,  Xa=R,6x,  u,.u,u,  v,  =  (u./Rj)v,  Wa=u„w,  Pa=p.p,  Pd=P«u» 
1,  The  index  "d"  corresponds  to  dimensional  variab¬ 
les,  index  -  to  the  flow  in  the  initial  section  at  r  - 

».  Ej=Po,u..6/p.is  the  Reynolds  number. The  analysis  of  the  Na- 
vier-stokes  equations  at  the  limit  r-O(l),  x  0(  ) 

lows  to  find  the  first-approximation  equations; 

5u  6u  6p  1  5  6u 

p  ( u —  +  V — )  =  "  ~  ^  ^  ' 

6x  6r  6x  r  6r  6r 
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5g  5g  5  i  6g  2  6/a 

p  (u —  +  V — )  =  r  (/I  )  -  g 

5x  6r  6r  r  6r  r  6r 


6(pu)  1  6(pvr) 


6x 


+  —  - - 

r  6r 


=  0 


1  6p  ’T 

_  —  =  p- 

r  6r 


6h 

6p 

6p 

1  6  6h 

V— ) 

=  U -  ^ 

-  V -  + 

- (-  r  — )  + 

6r 

6x 

6r 

r  6r  6  6r 

-  6u 

1  6g  g 

+  M 

{—f 

+  ( - -  2— )2 

6r 

r  6r 

(1.4.1) 


^  P 

h  = - ,  M  =  M(h)  . 

K-l  p 

where  g=wr . 

The  equation  for  h  may  be  changed  by  the  equation  for  p: 


6p 

Sp 

6u  1 

6(vr)  1 

u —  + 

V -  = 

-  ')(P(—  +  — 

- )  +  (K-D- 

6x 

6r 

6x  r 

6r  r 

6r  6  6r 


(1.4.2) 


r  6u 


+  (K-i)M 


1  6g 


(_)2  ^  ( - -  2— )‘ 


6r 


r  6r 
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The  system  (1.4.1) 
mulate  the  problem 

potential  flow. 


is  a  nonlinear  parabolic  system.  We  for- 
statement  for  the  vortex  surrounded  by 


initial  profiles  ,  g.  .  h,  are  deterBined  at  the 

Uion  x=0  (u,r,=l,  The  pres¬ 
sure  profile  Pi(r)  is  restored  by  Integration  with  the  boun¬ 
dary  condition  ■  If  the  external  pressure  dist¬ 
ribution  p(x,»).p.(x)  is  deternined,  the  axial  velocity  and 
enthalpy  at  infinity  are  found: 


h(x,oo)  =  he  (X)  - 


(K-DM. 


Pe 


(K-D/K 


(1.4.3) 


(x,co)  =  Ue^  (x)=  1  + 


(K-l)M. 


(^-D/'y 


l-(1(M^  Pe  ) 


Then 


the  boundary  conditions  are  formulated  as  follows 


u(0,r)  =  Ui(r),  g(0,r)  =  gi(r).  h(0,r)  -  h^ (r) 
u(x,w)  =  Ue(x),  g(x,oo)  =  g^,  h(x,co)  =  he(x); 


(1.4.4) 


g(x 


,0)  =  0,  lu(x,0)|  <  00,  lh(x,0)l  <  00,  v(x,0)  =  0  . 


An  interesting  example  of  the  problem  is  the  development  of 
the  potential  vortex  (1.3.9)  caused  by  viscous  effects.  The 
choise  5  =  r.  results  in  the  following  boundary  conditions: 


Ui  (r)  =  1,  gi  (r)  =  gw 


M. 


2  n 


L^-l 


1/2 
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hi 


=  h^(l  -  — ) .  Pi  (J^)  =  Pc^  ^  " 


1 


,  r  >  1,  (1.4.5) 


(x)=l/  he(x)=const,  Pg(x)-const 


problem  at  fixed  gas  model  (K.e,/!)  contains  the  only  pa¬ 
rameter  -  However,  this  problem  has  a  singular  inner  bo 

undary  -  thi  vacuum  core,  and  requires  special  treatment. 

The  equations  (1.4.1)  may  be  represented  in  the  form  resol¬ 
ved  relatively  x  -  derivatives.  The  procedure  is  similar  to 
those  used  in  derivation  of  (1.3.5),  (1.3.7)  results  in. 


f  +  ( - )  '  + 


y  Kp(i-m^) 


(1.4.6) 


1  5u 


=  -  f - (f—  +  f  '  )/(l-M^  )  + 


Dg/u  -  Di 


Kp(i-m^ ) 


1  5p 


p  6x 


=  _  f - f  -  +  (f—  +  f  '  )/(l-M^  ) 


Dg/u  -  Di 


Kp(i-m^ ) 


1  5p 


—  =  (f—  +  f  )/(i-mM 


pu^  6x 


Di/M^  -  D3/U 


'Kp(i-m^  ) 


(1.4.7) 


1  5g 


g'  D2 


g  5x 


g  ug 
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6  6u  1  j-  6  6g  5M 

=  — (2My— )  ,  D2  =  -  2y— (M— )  -  2—  g 
5y  5y  p  L  6y  6y  6y 

6  fi  6h  r  6u  6g  g 

D3  =  —(2-  y— )  +  (K-i)M  2yi—f  +  (—--) 
6y  6  6y  '■  5y  6y  y 


(D3/u)-Di 


p'  rDg/u-Dn  pgDg 


^  1-M 


o 

2y  u 


where  f=(vr)/u,  y-r^ /2  and  the  primes  denote  the  y-derivati- 
ves.  The  boundary  conditions  for  the  first  equation  are: 


f (X, 0 ) =0 , 


f-(x,c«)  =  (1-Me^) - /(Pe^e) 


If  the  profiles  u,  g,  h,  p  are  known  at  some  section,  the 
coefficients  of  (1.4.6)  are  determined.  The  problem  for 
(1.4.6)  is  inhomogeneous  and  coinside  at  f  ('^)  0 ,  Dq  0  wi 

the  problem  (1.3.7)  at  a=0.  Therefore,  if  the  spectrum  of 
the  problem  (1.3.7)  does  not  contain  a=0,  the  problem 
(1.4.6)  has  unique  solution,  the  right-hand  sides  of  (1.4.7) 
can  be  calculated  and  the  profiles  may  be  continued  on.  In 
opposite  case,  the  problem  has  singularity  or  bifurcation  at 
the  section  and  the  solution  either  can  not  be  continued  on 
at  all,  or  may  be  continued  in  various  ways. 

At  presence  of  sonic  point  M=l,  the  special  correlations 
must  be  satisfied  in  it  to  avoid  the  singularity  in  the  dif¬ 
fusion  terms.  In  this  case,  the  equations  (1.4.6),  (1.4.7) 

may  be  found  unsuitable  for  numerical  implementation,  and 
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original  statement  (1.4.1),  ( 1  ._4 . 4  )  with  implicit  method 
of  integration  is  to  be  used. 

The  equations  (1.4.7)  may  be  also  used  for  description  of 
the  inviscid  vortex  development  under  the  action  of  the  ex¬ 
ternal  pressure  changes.  In  this  case,  the  axial  scaling  is 
determined  by  the  scale  of  the  external  pressure  pg  and  all 
diffusion  terms  are  to  be  discarded.  This  approach  is  equi¬ 
valent  to  the  solving  of  the  nonlinear  problem  (1.2.7) 
step-by-step  from  the  initial  section. 


2.  The  structure  of  flow  at  the  initial 
region  of  vortex 

Expressions  (1.4.5)  cannot  be  used  directly  as  initial  con¬ 
ditions  for  numerical  integration  of  equations  (1.4.1.)  be¬ 
cause  of  presence  of  the  singular  vacuum  core  at  0  <  r  <  1. 
Near  it's  boundary  relation  pu/fX  decreases  so  that  effects 
of  viscous  forces  are  significant  at  small  distance  from  the 
initial  section.  Hence  singular  layer  appears  near  the  cy¬ 
lindrical  surface  at  r=l .  This  layer  reveals  essential  pro¬ 
perties  of  the  classical  mixing  layer  but  in  our  case  it  is 
more  properly  to  name  it  "penetration  layer" . 


2.1.  Representation  of  the  solution 
in  the  penetration  layer 

Let  us  consider  the  flow  in  the  neighbourhood  of  the  vacuum 
core  and  introduce  new  variable: 

r^  -1 


(2.1.1) 
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Of  the  solutions  are  represented  as 


u  =  1  +  + 


g  =  g^+  91 


+  .  . 


h  =  hi  (U  + 


p  =  XPi  (U  + 


•  •  ♦  / 


(2.1.2.  ) 


p  =  x^^K  Pi  (U  +  • 
vr=  x~  ^  ^  'i  Vi  (^  )  + 


M 


=  h:i  hi  + 


where  the  dots  denote  higher  terms. 

substituting  (2.1.2)  in  the  equations  (1.4.1)  and  performing 
the  limit  x  0,  we  obtain 


pl  [  -  (Ui -^1  u' )+Vi  u/  ]  = 

1 


^-1 


-Pl  +  -  iPi  '  +  (hiUi 


^-1 


Pl  [• 


(gi-Ui')  +  Vigi']  =  h-^((higi')'  -  2hi'gi), 
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K-l 

Pi  [ -  (hi  -  ^hi  '  )  +  Vjhi  ’  ]=  Pi 

1 


K-i 

—  ^Pl  '  +  ViPi  '  + 

K 


1 

—  h;i(hihi-)', 
6 


1  1-1 

—  pi - ^Pi  '  +  (PiVi)'  =  0, 

1  1 

1 

Pi  '=Pig^  '  Pi^i  =  Pi  ^ 

1-1 


(2.1.3) 


where  the  prime  denotes  derivatives  on  variable 

System  (2.1.3)  is  the  system  of  ordinary  differential  equa¬ 
tions  with  boundary  conditions  obtained  from  matching  with 
the  functions  in  the  exterior,  i  +*«,  and  interior,  ^ 
fields.  We  suppose  that  Pi(i)  tends  to  zero  as  ^  -<=«  more 

quickly  than  1/  i  ,  i-e.  pi^  ^  0,  ^  Than,  from  the 

latter  three  equations  (2.1.3)  we  obtain 

i  1 

PlVi=Pi  (—  -  -  )  .  (2.1.4) 

hi  gfo 


Using  (2.1.4)  we  see  that  (2.1.3)  may  be  written  as 


6 

(hihi  '  )  '  +  -  Pi  (hi  '-g^)  =  0, 

2 
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K 

hiPi '  -  —  Pi  =0^ 

1-1 

gL 

(hi  Uj  '  )  '  +  hi  Ui  '  -  Pi  ( 

2  hi 


(2.1.5) 


Ui  ' 


-  gi  — ) 

hi 


0, 


gi,  gi 

(hi  gi  *  )  hi  '  (gi  ~25go)  Pi  ( 

2  hi 


1 


gi  ’) 


«2 


0  . 


The  solution  of  the  first  and  second  equations  in  (2.1.5) 
is  found  by  means  of  confluent  hypergeometric  functions  M,  U 
[42]  : 


6  1-1  1  2^-1  6  ')(  -1  Pi 

hi  =  exp[- - Pi  ]  (CiM(  ,  1,  )  + 

2  gi  1  2  •)(  g^ 

21-1  6  1-1  Pi 

+  C2U(  -  ,1, - )),  (2.1.6) 

1  2  ')(  gi 

were  Ci  and  Cg  are  arbitrary  constants. 

We  consider  that  Pi  0  at  ^  Therefore,  for  small  va¬ 

lues  of  Pi  from  (2.1.6.)  we  obtain 

6  1-1  1 

hi  =  q  (1+ -  ( - )^  -  Pi  +  0(pi^  ) )  - 

2  1  g^ 
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C2  6  y  -1  1 

- - [ln( - - Pi  )  +  0(1)],  (2.1.7) 

r  ((2](-l)/K)  2  K 

were  F  is  the  Gamma  function. 

If  Cg  ^  0,  expression  (2.1.7)  contains  logarithmic  singula¬ 
rity: 

-Cg  6  K  -1  1 

hi  ~  ln(  Pi  )  +  ...  , 

r((2')(-l)]()  2  ](  gt. 

and  hence  the  asymptotic  behavior  of  the  solutions  of  first 
and  second  equations  in  (2.1.5)  at  ^  is 

21  K  -1  i-  io 

Pi  ~ - gLexp[g^(2r  ( - )  -  )  ], 

6  1-1  1  Cg 

(2.1.8) 

C2  1-1  i-io 

hi  ~  -  g«(2r( — ) - )  ^ 

r((2i-i)/i)  y  cg 

where  is  the  constant  of  integration.  From  (2.1.8)  we 
find  that  hi  ^  +0^  at  ^  so  that  heat  flux 

1  5h  1  1^1 

M - -  —  hi  hi  '  ~  2C2  ( — )  - 

r  6r  h  1-1  r((l-l)/')() 


tends  to  the  constant  value. 

The  absence  of  the  heat  flux  inside  the  vacuum  core  means 
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that-i”  (2.1.6.)  it  is  necessary  to  assume  €2=0.  In  this  ca¬ 
ge  at  ^  the  main  terms  in  asymptotic  expansions  for  hj  , 


6  K-l  ^  1 

hi  ~  Cj  ( 1  +  — (  )  Pi  )  , 

2  Y  g - 

(2.1.9) 

If  i-io 

Pj  ~  exp[ -  g*  - ]  •. 

If-l  Cl 


Let  us  investigate  the  behavior  of  the  function  gi  (O  at  ^ 
-00.  After  considerable  manipulation  we  obtain  from  (2.1.5) 
the  equation 

d^gi  1  dgi  1  If-l  2  K-1  dhi 

Pi  +  (1  +  Pi  )  -  (  )^9i  -  ~  • 

dpi^  2gL  K  dpi  2gi  K  g«,  “)(  dpi 

(2.1.10) 

The  general  solution  of  (2.1.10)  is  represented  in  the  form: 

00  00 

gi  =  I  Akp/  +  Inpi  I  Bj^p/  ,  (2.1.11) 

k  =  0  k  =  0 


Were  are  constants  to  be  found. 

If  Bq=o  in  (2.1.11),  the  solution  .of  the  equation  (2.1.10) 
tends  to  constant  at  Pi  0 : 

26  K-1 


g®?  ”)( 


gi 
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asymptotic  behaviSr  of  circulation  does  not  match  with 
boundary  condition  at  r=0 . 

The  second  solution  of  (2.1.10)  has  a  logarithmic  singulari¬ 
ty 


gi  ~  Bolnpi  , 

that,  using  (2.1.9),  reduces  to  linear  dependence 

gi  ^  Bog^ - .  (2.1.12) 

^-1  q 

Substituting  (2.1.12)  into  the  expansion  (2.1.1),  (2.1.2)  we 

find 

r^  -1 

g(x,r)  =  g^+  91  (x^  ^  ^K)  -  )  ~ 

2 

")(  r^  -1 

^  gu+  Bog2^ - . 

1-1  2q 


Boundary  condition  g(x,0)  =  0  determines  the  coefficient  Bq : 

2q  ')(-i 

Bo  - - . 

1 

Therefore,  the  circulation  distribution  in  the  initial  sec¬ 
tion  of  the  vacuum  core  is  following 

,  at  0  <  r  <  1  (2.1.13) 

Indeed,  such  circulation  means  that  the  vacuum  core  is  spin- 
^^ng  around  the  vortex  axis  like  a  solid  body.  Its  angular 
^®locity  of  rotation  is 
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6wd  1  2 

—  =  -  (  — — 

5rd  Ko  'K-1  ^d* 

Function  Ui(^)  satisfies  the  equation  that  may  be  obtained 
from  (2.1.5): 

Uj  1  ^ 

pj -  +  (1  + - Pi  ) -  -  - ( - = 

dp/  2g^  ")(  dpi  2g^  'K 

1  ^-1 

=  - ( - )^(hi-gj.  i)  ■  (2.1.14) 

2g^  1 

Left  hand  —  side  in  (2.1.14)  is  similar  to  that  in  (2.1.10). 
The  solution  of  (2.1.14)  at  is  expanded  as  a  series: 

00  00 

Ui  =  I  Dkp/  +  inpi  I  E^p/  (2.1.15) 

k  =  0  k  =  0 


Absence  of  an  infinite  term  in  (2.1.15),  Eg  -  0,  means  that 
in  the  whole  region  0  C  r  C  1  the  axial  velocity  slightly 
differs  from  the  initial  value.  The  difference  between  them 
tends  to  zero  as  x  0 :  u(0,r)  -  u(x,r)  =  0(x^')(  ^  ”)( )  •  In 
this  case  the  total  enthalpy  hQ(0,r)  in  the  initial  section 
of  core  is  a  parabolic  function  of  r: 

ho(0,r)  =  1/2  +  h^  r^  ,  at  0  <  r  <  1.  (2.1.16) 


In  contrast  to  this,  if  the  vortex  generator  does  not  produ 
ce  considerable  changes  of  the  total  enthalpy,  we  must  take 
into  account  the  initial  axial  velocity  profiles  different 
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onstant  At  constant  initial  total  enthalpy  pcofile 
from  ^  circulation  distribution  obtained 

-f  :3rt./::ia::eLcUv 


nu  r  ^  \  ^  ^  ,  0<r<l- 

u(0,r)  =  (!■*■  2hgo( 

Hatching  (2.1.17)  with  (2.1.15)  yields 


(2.1.17) 


(2.1.18) 


Eo  =  " 


The  unknown  constants  q  ,  U  (2. 1-9)7 

ace  determined  from  matching  the  solutions  o  .1-)  -th 

r::rM" ^r;a::ior:/:;::gic™etric  tunction 

purpose  w  , ,  ov  2  x  x  n  oo  [  42 1  .  Then,  from 

N((2’X-l)/-«,  1.  6(T(-1)  Pl/<2Tl9>e>>  =>’'  P‘ 

(2.1.6)  we  find: 


c,  6  ^-1  1 

- [ - —  Pi  ] 

r((2K-i)/^)  2  K  gL 


(2.1.19) 


substituting  (2.1.19)  into  the  second  equation  (2.1.5), 

obtain 


2  1 

Pi  ~ - 

6  “)(  -1 


g2  2^-1 

!_r( - )  (^^l)  ,  (2-1-20) 

c,  t 


where  L  is  arbitrary  constant. 

D  ^  /0111  f2l4)  (2.1.19)  and  (2.1.20)  cons 

Erom  (1.4.5),  (2.1.1)7  (2.i.4j,  i 

tants  Cl  and  are  determined: 
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Cl 


K-1 

h„,r( - ) 

1 


-  K-1- 

-  4  - 

Y  - 

-  6  - 

(2.1.21) 


a  result,  at  ^  “  we  have: 


~  2hatf^ 


Pi 


^_1  l/(^-l) 

gi — 2  K 

1 


1  /  (  v-  1  )  1  /  (  V”  ^  ^ 


Pi 


(2.1.22) 


2.2  Numerical  study  of  the  flow  at  the  initial  region 

The  finite  difference  method  was  used  for  integration  equa¬ 
tions  in  (2.1.5)  to  calculate  vortex  development  at  the  ini- 
--tial  region  of  flow.  The  independent  variable  i  was  conside- 
-red  in  the  range  -3.5  <  ^  <  3 . 5 .  The  second  order  method  of 

^“Keller  was  used  to  approximate  derivatives.  At  the  first 
;  stage  the  Caushe  problem  for  the  first  and  second  equations 
r-  in  (2.1.5)  was  solved  with  initial  conditions  (2.1.9)  stated 
;  at  ^  =  -3.5.  The  form  of  equations  (2.1.5)  does  not  change 
if  an  arbitrary  constant  is  added  to  independent  variable  . 
Therefore,  constant  (2.1.9)  may  be  chosen  so  that  = 

0  in  (2.1.20).  On  each  step  of  the  iteration  procedure  cons¬ 
tant  was  being  found  from  (2.1.20)  by  solving  the  Coushe 
problem  with  given  value  .  Then  the  value  of  chan 

ged  to  provide  the  condition  =  0.  The  process  was  finis¬ 
hed  when  appropriate  accuracy  was  reached.  It  is  clear  that 
such  iteration  procedure  is  necessary  because  of  desagree- 
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f  between  asymptotic  character  of  expansions 
’  1,22)  and  the  finite  interval  for  independent  variable. 

,^on  as  functionspi  (U  and  h,  (U  are  found  the  latter 
"  eduations  in  (2.1.5,  are  solved.  These  equations  are  ll- 
with  respect  to  u,  and  g.  ,  and  no  iteration  procedure 
required  to  satisfy  the  boundary  conditions.  On  the  rig  ^ 

*  „ndary  of  interval  conditions  Uj  =  0,  gj  -  0  were  consl  e 
a  on  the  left  boundary  the  derivative  g'  was  given:  g,  eg„ 
■'^unction  u(U  was  calculated  in  two  cases  corresponding  to 
two  boundary  conditions  on  the  left  boundary  of  the  interva 


Uj  '  =  0 


(2.2.1) 


(2.2.2) 


The  results  of  the  calculations  of  functions  h,  (U  -  P.  (U  , 
p.(U,  9l(t)  at  6  =  0.7.  H„=  3,  K  e  1-4  are  presented  in  Fl- 

gures  42-45  respectively* 

solution  U,,(U  satisfying  the  condition  (2.2.1)  is  repre¬ 
sented  in  Fig. 46.  The  solution  u.2(U  for  condition  (2.2.2) 

is  shown  in  Fig. 47. 

Representations  (2.1.1),  (2.1.2)  and  functions  (i)  ,  Pi(Ut 

91  (U,  ui  (O  obtained  above  determine  behavior  of  the  flow 
in  thi  initial  region  of  vortex.  Thus,  pressure,  density  and 
enthalpy  increase  monotonously  inside  the  vortex  as  axial 

coordinate  x  increases. 

The  enthalpy  on  the  vortex  axis  has  non-zero  value: 
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h(x,o)  ~  - 

K-1  K-1 


( — ) 
6 


.  (2.2.3) 


Returning  to  dimensional  variables  we 
the  following  correlation  of  vortex  parameters: 


obtain  from  (2.2.3) 


hd(Xd/(R5S),0) 


Xd 


( 1  -  v)  / 


V  t 


C7*» 


4  - 

(-y- !)/')(  1  “X-l 

-K-i- 

-  r( - ) 

*-61^- 

'i  ^ 

(2.2.4) 


particularly,  at  ^  =  1 • 4 ,  6  =  0 . 7 ,  (2.2.4)  reduces  to 

hd(Xd/(R  6),0)  7.56 


(Xd  /RgS) 


0.29 


M. 


in  the  case  when  the  Initial  axial  velocity  profile  is  con 
tant  condition  (2.2.1)  results  in  negative  values  - 
tion  u,,(U  (eee  Fig. 46).  Therefore,  the  flow  is  retarded 
least  for  that  distance  from  the  initial  section  where  ex 

pension  (2.1.2)  is  valid. 

If  the  initial  total  enthalpy  profile  is  constant,  i^e.  the 
initial  velocity  distribution  is  given  by  (2.1.7),  e 
tardation  of  the  flow  is  less  significant  and  velocity  on 
the  axis  stays  constant  in  the  first  approximation: 


u(x,0)  =  (1  +  2h^) 


1/2 
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prom  expansions  (2.1.2)  and 
ui  all  initial  profiles 
cedure  may  be  found  . 


calculated  functions  ,  Pi  , 
required  for  the  numerical  pro- 


3. 


numerical  solution  of  the  quasi-cylindrioal 
approximation  system  for  compressible  flow 


(1  4  1)  is  non-linear  parabolic,  just  as  the  system 
Tboundary  layer  equations.  Effective  numerical  methods  to 
solve  such  equations  with  given  initial  and  boundary  con  i 
tions  are  known.  All  programming  and  computing  was  per  orme 

using  VAX-780  computer. 


3.1.  Problem  statement.  Numerical  method 

The  system  (1.4.1)  was  solved  numerically  using  verlables^x. 
y  =  r^/2  inside  the  region  x  >  0, 


0  <  r 


The  upper 

boundary  Condition  was  stated  at  point  r™,,  =  10.  To  avoid 
the  errors  caused  by  finite  region  of  calculation  the  folio- 
wing  asymptotic  expans ions were  used: 


h(0,rmax )  =  (I' 


)> 


^  max 


P(0fr-max  )  =  Pe  (1" 


„2 

^  max 


(3.1.1) 


u(0,rmax)  = 


g(  0  ,rniax  )  ■ 


solution 


must  be  bounded  on  the  axis 


and  it  results  in 


jnditio’^s  [43]: 


^  0, 


^  0, 


r  0 . 


(3.1.2) 


M  4  1)  with  conditions  (3.1.2)  the  fol- 

l/.un.ucns  u  ana  .  on  a.is  .s  .onna: 


1  6h 


[p  h 


-//(/-  1 )  1  + 


//  (7'- 1 ) 


=  0, 


(3.1.3) 


d  u 


dx  2 


+  h  ] 


2JLI  [- 


1  6h 


6r  6r 


1  6u 


=  0. 


in  addition  to  the  conditions  (3.1.1) 
ons  for  V  and  y  on  the  axis 

v(x,0)  =  g(x,0)  =  0 


(3.1.3)  the  conditi- 


(3.1.4) 


are  satisfied. 


I  r\  Y-  \  afOr)  must  be  spe 

Initial  profiles  u(0,r),  h(0,r  ,  P  j  consider  two  types 


nroblem  statement.  we  consxi.t=.  - - 

cified  for  correct  pr  different  condi 

iifO  r)  corresponding  to  aiiietei 
of  initial  profiles  u(  , 

tions  of  vortex  generation 


u(0,r)  =  1 


0  <  r  < 


(3.1.5) 
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and 


u(0-,r)  =  1,  r  >  1  , 

u(0,r)  =  [1  +  2h„(l-r2 )]*'",  0  <  r  < 

initial  profiles  of  h,  p,  g  are  the  same 

sidered: 


1  .  (3.1.6) 

in  both  cases  con- 


p(0,r)  =  h(o,r)  =  0, 
g(0,r)  =  g^r^ 
at  0  <  r  < 

P(0,r)  =  p^(l 

h(0,r)  =  h^(l 

g(0,r)  =  g^ 
at  r  >  1 . 

The  absolutely  stable  finite-difference  method  with  se¬ 
cond-order  approximation  of  the  radial  derivatives  and  vari¬ 
able  grid  size  in  radial  direction  was  used  for  numerical 
solution  of  equations  (1.4.1).  A  second  order  implicit  sche¬ 
me  was  used  for  approximation  in  axial  direction. 


(3.1.7) 

1 

r^ 

1 

— ). 


1' 

! 


ill 


3.2  Numerical  results 

■  tial  profiles  were  found  usin?  asymptotic  expansions 
"  ,712,  The  profiles  of  functrons 

^Inrx  .rvJ,  9(x,.r)=g(x,.r,/9.,  P(x,  .r)=P(x,  ,r,/p 

.  .nown  in  figure  52  at  Pne  same  vai_  o 

:u:tro:tiorof  axiai  velocity  m  tne 

:::  of ''Constant  initial  total  enthalpy  ie  demonstrated  rn 
■igure  53. 

4-  rn-F  the  nrofiles  downstream  from 
me  subsequent  development  .  ^  _3  Y-i  4  6=0.7  and 

•  section  was  calculated  at  M^-3,  K-1.4,  u 

the  initial  section  w  viscosity  |1  was  consi- 

— =:»  rs. 

r  rr,,.!.".... -  ""t:.  »:.tr 

r-::rjnrr;r.ru:“.rr.::;.. ...... 

at  Reynolds  number  -  2  10  . 


2.510  ^ 

5.5-10  ^ 

9.5- 10"^ 

CO 

O 

CM 

1  ^ 

1.110'^ 

2.0-  10^ 

0 . 29 


6.0-  10- 


Xd  /O 


8.9-  10" 


1.5-  10^ 


2 . 4  -  lO'- 


The  profiles  h,  p,  p,  'g ,  u,  hg  at  various  axial  positions 
are  demonstrated  in  figures  54-59.  The  profile  of  the  total 
pressure  Pg  calculated  taking  into  account  the  local  Mach 
number  is  shown  in  figure  60.  All  these  functions  are  scaled 
by  their  values  in  the  external  flow  at  r-^oo.  Radial  and  azi¬ 
muthal  components  of  velocity  (v,w)  and  local  Mach  number  M 
ate  represented  in  figures  61-63.  Swirl  angle  (p=arctg(w/u) 
between  velocity  vector  and  the  axis  of  the  vortex  is  shown 
in  figure  64  ((p  is  measured  in  degrees)  . 

The  numerical  results  represented  in  figures  54-64  were  ob¬ 
tained  at  initial  conditions  correspondent  to  constant  axial 
Velocity  at  the  initial  section.  The  results  demonstrate 
^bat  the  axial  velocity  inside  the  vortex  is  reduced  near 
the  initial  section.  This  tendency  is  obvious  in  figure  65 
(curve  1),  where  the  axial  velocity  on  the  axis  is  shown. 
The  retardation  of  the  axial  flow  is  caused  by  increase  of 
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.sure  downstream  from  the  initial  section.  The  reason  of 
pressure  increase  is  diffusion  of  gas  having  small  azi- 
^^thal  velocity  towards  the  vortex  core  resulted  in  reducing 
the  radial  pressure  overfall.  The  distribution  of  the  to¬ 
tal  pressure  on  the  axis  is  represented  in  figure  67. 

numerical  calculations  were  performed  also- for  the  case  of 
onstant  total  enthalpy  in  the  initial  section.  It  was  found 
that  the  distributions  of  h,  p,  g  are  quite  close  to  those 
obtained  at  constant  axial  velocity  in  the  initial  section. 
Significant  difference  is  observed  only  for  axial  velocity 
and  total  enthalpy  on  the  relatively  small  distance  down¬ 
stream  of  the  initial  position.  The  distributions  of  axial 
velocity  u(x,0),  scaled  total  enthalpy  ho(x,0)  and  total 
pressure  Ps(x,0)^  Mach  number  M(x,0)  along  the  vortex  axis 
are  shown  in  figures  65-68.  The  curves  1  correspond  to  the 
case  of  constant  axial  velocity  in  the  initial  section  and 
curves  2  -  to  the  case  of  constant  total  enthalpy.  It  could 
be  seen  that  distributions  of  the  total  pressure  in  both  ca¬ 
ses  almost  coinside.  The  retardation  of  the  flow  inside  the 
vortex  found  in  the  first  case  was  observed  in  the  second 
case  too.  However,  the  difference  between  the  flows  caused 
by  different  initial  conditions  decreases  quickly  downstream 
-from  the  initial  section.  It  is  interesting  that  in  the  se¬ 
cond  case  the  total  enthalpy  quickly  drops  almost  to  the  va¬ 
lues  found  in  the  first  case.  Moreover,  positions  of  minimal 
axial  velocity  and  total  enthalpy  on  the  axis  are  rather 
close  in  both  cases. 

It  may  be  concluded  that  development  of  the  vortex  downstre¬ 
am  from  the  initial  position  is  caused  by  two  diffusion  pro¬ 
cesses.  The  first  process  is  spreading  of  the  gas  towards 
the  vortex  core.  This  process  determines  the  profiles  at  re¬ 
latively  small  distance  from  the  initial  section.  The  second 
process  is  regrouping  of  the  profiles  under  action  of  visco¬ 
sity. This  process  is  relatively  slow  and  determines  profiles 
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large  distance  from  the  initial  section.  Two  functions 
h  Characterize  radial  size  of  the  vortex  were  conside- 
The  size  5„  was  determined  as  the  radial  position  where 
velocity  u=0.99.  The  second  size  6,^  was  detemined  as 
the^radial  position  where  the  swirl  angle  tp  is  maximal  (fi¬ 
gure  64)  . 

nistributions  of  these  values  along  the  vortex  is  demonstra¬ 
ted  in  figure  69.  It  could  be  concluded  that  the  vortex  de¬ 
velops  smoothly  and  its  radial  size  monotonously  increases. 


4.  Solutions  describing  the  flow  far 
downstream  from  initial  section 

The  asymptotic  expansion  for  dissipating  incompressible  tra¬ 
iling  vortex  far  from  aircraft  was  obtained  in  Batchelor's 
work  [29].  It  was  shown  by  Batchelor  that  an  arbitrary  cons¬ 
tant  appears  in  his  asymptotic  expansion.  We  extended  the 
expansions  and  found  terms  which  contain  two  arbitrary  cons¬ 
tants  else. 

4.1.  Generalization  of  the  Batchelor's  expansion 
Let  us  introduce  a  new  independent  variable: 


The  asymptotic  expansions  of  the  solution  are 
Inx  1 

U  =  1  +  -  Uj  (T\)  +  “  ^2  ' 

X  X 


^  1  -  exp(-Ti)  .  —  gi(n)  -  ' 


=  1  + 


hi  (n)  +  “  ^2 


=  1  + 


(n)  +  "  P2 


(4.1.2) 


=  1  + 


pi  (K\)  +  "  p2 


h  +  -  Vg (n)  +  • • •  ' 


M  1 


(- -  hi  +  "  hg  +  •  •  •  )  ' 


dh  1  h=hoo  ^ 


where  the  dots  denote  terms  which  have 


the  form  0(x*^(lnx)  ) 


for  all  k ,  1  • 


,  wv  (4.1.2)  into  equations  (1.4.. 
After  substitution  (  •  •  matching  we  find 
in  the  first  approximation  as  a  resul 


-  —  a^..  exp  (  -11) 
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1 

=  -  -  1f\  exp(-lfl)  [g^6jloc(Gi  2  (0)  "^1 1  ^ ' 

8 


1 

=  -  6  exp(-6lf\) 
4 


(4.1.3) 


1 

=  0  ,  Pi  =  -  “  6  exp(-6n)  . 

4 


Vi 


1  1 
-  —  1f\[6  exp(-6T\)  +  “  gL  exp(-T\)]  ^ 


where 


jii 


dfl 


dh 


h— hoo 


and 


Gii (n) 


n 


d1f\ 


/exp(-6n)  [exp(1fl)-n-l 

0 


Gi2(n)  = /ncii  (TDexpc-iDdn 

n 


Arbitrary  constants  Cy ,  Cg , 


Cj^  appear  in  the  second  terns . 
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U2 


92  “ 


h2  = 


P2  = 


P2  - 


V2  = 


1 

-  exp(-6t\)  [UCIfl)  +  Cu  ]  , 

8 

1 

-  'r\  exp(-1f\)  [6)iaj(gdG2 1  CU)  +  4CjjGii(Tl))  + 
8 


+  (g^/2  +  46/(1  +  6))G2  2(11)  +  1 


(4.1.4) 


1 

-  6  exp(-6T\)  [11(11)  +  Ch  ]  ' 

2 


1 

- p(iri)  , 

4  K-i 

1  1 

- P(Tf\)  -  -  6  exp( -6tI)  [H(T\)  +  2Ch/g£]  , 

4^-1  2 


1  1 
-  g^  [l-exp(-'ri)  ]  +  -  [l-exp(-6n)]  + 


+  “  g^ll  exp( -1^)  [U(T[)  +  Cu  ]  - 

4 


1  1 

-  6n  exp(-6Tl)[H(n)  +  2Cji/g^]  - - T\P(n)  ^ 

2  ^-1 
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Where 

00  d’fl 

P(n)  =  /[I  -  exp(-ifi)  f  —  > 

T\ 


U(Tl)  =  /CIIPCA)"!  +  exp(-T\)  )exp(Tl) 


1  n 

H(n)  =  “  /  [exp(-2ifi)  -  exp(-6n)  "  np(Ti)  + 

2  0 

dTl 

+  2(1  -  exp( -tl)  /ll]exp(6t\)  , 

n 

n 

Ggicn)  =  /[2H(n)(exp(n)-  i  -  n)exp(-6n)  + 

0 

dn 

+  gi2(o)(i  +  n)  -  Gi2 (n)exp(-n) ] ]—  , 

^2 

n 

*^2  2  (^)  “  ^  ^  ] 

0  yf 

Using  the  exponential  integral  Ei(n)  [42]  we  may  produce 
nie  calculations: 

P(0)  =  2ln2, 


so- 


Gj2(o)  =  6/(i+6)+(i-6)in( (i+6)/6)  , 
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to  obtain  correlations 


h(>^ 


0 )  ~  hoo(  ^ 


6  Inx  6 

_  -  +  C{j  ) 


4  X 


2X 


9<x? 


u(x,0)  ~  1 


_  -  (  Inx  -  Cy  )  / 


8X 


p  (  X  ,  0  )  ~  Po>  ( ^ 


t  l’^2 


K-1  2X 


(4.1.5) 


1  6g(x,r) 


6r 


Inx 


■[1 


-=0  2X 


8X 


1  +  6 


,.6  gi 

+  (1-6)  In-  +  + 

6  26  1-6 


) )  +  Cg  ]  . 

8X 


4.2  comparison  of  the  numerical  results  with 
asymptotic  solutions 

with  asymptotic  solutions  it  is  neces- 
erform  comparison  with  asy  P  ^  Correlations  (4.1.5) 

to  determine  used^or  this  purpose  so  as 

arge  enough  values  o  .h(x,0)  and 

and  asymptotic  solutions  of  u(x,0), 

(,0)/(r6r)  were  coinside  at  x=10. 

.  The  profiles  of 

alts  are  demonstrated  in 

ction  h.h/h„,  ^  ®  salid  curves  correspond  to 

73  at  various  axiai  h 
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numerical- profiles,  dashed  curves  correspond  to  the  first 
Approximation  (4.1.3)  and  dashed  curves  with  symbols  corres¬ 
pond  to  expansions  including  both  first  (4.1.3)  and  second 
(4.1.4)  terms. 

we  remind  that  the  first  approximation  does  not  depend  on 
constants  Cy  ,  Cj^  ,  Cg  . 


Distributions  of  axial  velocity  u(x,0),  scaled  enthalpy 
n(x,0)  and  pressure  p(x,0)  along  the  vortex  axis  are  repre¬ 
sented  in  figures  74-76  respectively. 

The  distributions  of  total  enthalpy  ho 0 )  number 
M(x,0)  along  the  vortex  axis  are  shown  in  figures  77-78, 
where  the  total  enthalpy  scaled  by  value  in  the  external 
flow.  Numerical  and  asymptotic  solutions  are  shown  by  solid 
and  dashed  lines  correspondently . 

A  good  agreement  between  numerical  and  asymptotic  results  is 
observed  at  large  values  of  x. 


5.  Comparison  of  the  numerical  results 
with  experimental  data 

The  comparison  was  performed  using  profiles  of  the  total 
pressure  and  stagnation  temperature  measured  experimentally 
in  two  radial  section  of  the  vortex.  At  the  first  step  nume¬ 
rical  and  experimental  profiles  were  correlated  in  the  down¬ 
stream  section.  The  ratio  of  the  total  pressure  on  the  axis 
(Y=ll3mm)  to  the.  total  pressure  on  the  conventional  boundary 
Pse  (Y=83mm)  was  found  from  the  Table  2,  and  then  the  nume¬ 
rical  data  from  figure  67  were  used  to  determine  position  x 
where  this  ratio  is  close  to  its  numerical  prediction.  The 
radius  of  the  vortex  was  approximately  determined  using  pho¬ 
tographs  taken  in  experiment.  Thorough  adjustment  of  parame- 


_  «as  performed  visually  qm  the  display  of  computer  to 
,  900d  correspondence  between  numerical  and  -  - 

'  fU«  TWO  parameters  were  adjusted:  reference  total 

ssure'and  initial  radius  of  the  vortex  vacuum  core  6.  Fi- 
'’"l  result  was  found  taken  into  account  possible  deflection 
"f/ne  Pbote  from  real  position  of  the  vortex  axis.  The  re- 
.It  is  represented  in  figure  79.  The  total  pressure  obtal- 
Ted  numerically  and  referenced  to  its  value  in  external  flow 
is  shown  by  solid  line  and  experimental  data  -  by  crosses. 
The  following  values  of  parameters  were  found:  inltia  radl 
„s  of  the  vacuum  core  6=5. 5mm,  the  ratio  of  the  total  pres 
s,re  on  the  vortex  boundary  to  the  total  pressure  in  exter- 
flow  Pse/P„=°.-3.  position  of  the  vortex  axis  V=114mm, 
and  position  of  the  section  considered  x  1.1  . 

Additional  information  about  the  flow  in  this  section  may  be 

obtained  from  figures  54-64  (curve  8).  In  particuar, 

dependence  of  the  swirl  angle  ,p  shown  in  figure  64  indicates 
that  maximum  of  the  angle  is  observed  at  radial  P°=3tion  ~ 
from  the  axis.  Analysis  of  the  experimental  data  (Tab  e 
results  in  rather  close  value  of  ^  and  radius  (ra=16mm)  whe¬ 
re  the  measured  total  pressure  Is  maximal. 

The  parameters  found  were  used  also  for  comparison  of  the 
numerical  and  experimental  profiles  of  the  stagnation  tempe¬ 
rature.  Both  profiles  shown  by  solid  line  and  crosses  cor- 
respondently  are  represented  in  figure  80.  The  stagnation 
temperature  decreases  toward  the  axis  and  quolltative  agree¬ 
ment  between  numerical  and  experimental  results  is 

Significant  discrepancy  is  probably  caused  by  relatively 
large  size  of  the  probe,  resulting  in  space  averaging  of  the 
experimental  data.  Furthermore,  both  profiles  are  rather 
close  to  unity  and  relative  discrepancy  is  small. 


Using  values  of  the  parameters 


found,  the  Reynolds  number  Re 


_2- 10^  and  dimensionless  coordinate  x . Correspondent  to  the 
first  section  where  the  profiles  were  measured  are  determi- 
However,  the  dimensionless  distance  Ax  between  the 
first  and  second  sections  turned  out  so  small  that  the  total 
pressure  practically  does  not  change  (Ax=4- 10"® )  .  On  the  ot- 
j^er  hand,  the  experimental  data  reveal  rather  appreciable 
change  of  the  total  pressure  between  two  sections  (Table 
i  2).  Therefore,  strong  disagreement  between  the  theoretical 
model  and  experimental  data  is  observed.  The  experimental 
data  indicate  that  radial  diffusion  is  really  more  intensive 
than  those  caused  by  molecular  viscosity. 

Two  reasons  may  be  considered  to  explane  this  disagreement: 

a:  Rather  high  level  of  acoustical  and  vortical  fluctuations 
is  observed  in  the  work-section  of  the  wind-tunnel.  The  dis¬ 
turbances  will  obviously  increase  the  diffusion  of  the  ave¬ 
raged  flow.  These  fluctuations  may  be  taken  into  account  by 
additional  "eddy  viscosity"  using  some  kind  of  turbulence 
model.  The  level  of  the  eddy  viscosity  will  be  evaluated  be¬ 
low  . 

b:  Additional  diffusion  may  be  caused  by  spiral  tangential 
discontinuities  in  the  vortex  which  were  not  taken  into  ac¬ 
count  within  axisymmetr ical  model  considered.  However,  radi 
al  distance  between  the  discontinuities  in  this  case  must  be 
very  small  (Ar<0.1mm).  Experimental  measurements  result  in 
spatial  averaging  of  the  profiles  on  much  bigger  distance. 
We  think,  the  discontinuities  can  not  completely  explain 
high  diffusion,  because  the  averaged  profiles  again  will 
fegroup  slowly. 

Diffusj^Qj-j  coefficient  was  estimated  by  comparison  of  numeri 
cal  and  experimental  profiles  of  the  total  pressure  in  two 
sections.  Numerical  and  experimental  values  of  the  total 
pressure  on  the  axis  coincide  if  the  first  section  is  situa- 
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l-ed  at  x=0.115. 

fjumerical  and  experimental  profiles  of  the  total  pressure  in 
the  first  section  are  shown  in  figure  81.  The  parameters  5, 
Pse/Po  were  not  adjusted  and  were  taken  the  same  as  in  figu¬ 
re  79.  Therefore,  only  value  x  was  specially  chosen  and  this 
is  equivalent  to  choise  of  the  diffusion  coefficient.  The 
agreement  is  not  so  good  as  in  the  second  section.  Indeed, 
the  first  section  is  situated  very  close  to  the  wing  and  the 
flow  can  not  be  consider  as  quasicyl indrical .  The  profiles 
of  the  stagnation  temperature  at  X=0.115  are  shown  in  figure 
82 . 


Taking  into  account  that  Ax=l . 1  -  0.115  ~  1,  effective  Rey¬ 
nolds  number  may  be  evaluated:  Re^  ^  ^  ~110mm/5 . 5mm=20..  Corres¬ 
pondent  diffusion  coefficient  /It  turned  out  very  high  in 
comparison  with  molecular  viscosity:  /Ij  ~  lO^/l. 
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Conclusions 


the  course  of  the  research  the  devices  and  models  neces¬ 
sary  for  vortex  investigation  and  its  interaction  with  i  - 
erent  bodies  were  developed.  The  vortex  investigation  and 
'  ^interaction  with  a  spherical  hlunt  body  at  Mach  number  3 

„d  Reynolds  numbers  from  1.1  to  3,7.10^  1/m  was  came 

out.  The  experimental  investigation  has  shown  the  following. 

1  prom  the  two  types  of  the  investigated  vortex  generators 
rectangular  wing-type  and  4-vanes-type  generators)  the 
nrst  type  of  the  generator  was  chosen  for  the  experiment^ 
me  rectangular  wing  produces  one  vortex  in  the 

part  Of  the  flow.  The  vane -type  generator  producee  ur^p^^ 

irs  of  vortices  (one  vortex  at  the  internal  * 

ends  of  each  vane) .  The  internal  v°btices  situated^clos 

each  other  flow  together  Into  one  vortex  while 

ternal  vortices  stay  apart  even  at  a  large  distance  from 

generator . 

3.  variation  of  the  vortex  diameter  with  the  ^ 

distance  from  the  generator  is  described  by  a 

,  .  ^4-  the  reference  vortex  diameter) 

on  d/dr.=f  (x/dp)  (Where  dr  is  the  reiere 

^  H  and  decree  of  the  flow  swirl  (the 

regardless  to  the  value  dp  and  d  g  n  , 

generator  angle  of  attach  «)  in  the  range  a  < 

tex  diameter  increases 

the  vortex  expansion  angle  is  abou 

le  increases  up  to  3  • 

3.  The  expanding  vortex  induces  in  a  supersonic  flow  a  weak 
compression  wave  as  well  as  an  expanding  viscous  wake  past 
the  wing  (vortex  generator). 


4.  The 


•  ^  at  =  20^  have  shown  that 

measurements  carried  out  at  a  n 
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.  -the  initial  cross  section  of  the  vortex  in  the  vicinity 
the  vortex  axis  a  region  of  the  low  stagnation  pressure 
°  is  formed.  The  minimum  value  Ps/Po  is  approximately  0.06 
’’f  the  free  stream  value.  At  a  distance  about  40d,  from  the 
°  itial  section  the  value  p./p,  at  the  vortex  axis  increases 
to  0.35.  The  laser  knife  pictures  also  show  the  fast  me- 
rease  of  density  in  the  vortex  core. 

,  The  stagnation  temperature  in  the  vortex  also  changes.  In 
the  initial  vortex  cross  section  the  stagnation  temperature 
,t  the  vortex  axis  is  lower  approximately  by  20  degrees  than 
in  the  free  stream.  However,  the  relative  value  of  the  stag¬ 
nation  temperature  Ts/T„  at  the  vortex  axis  in  the  initial 
section  is  about  0.90.  A  still  weaker  influence  of  the  flow 
swirl  on  the  relative  stagnation  temperature  is  observed  in 
the  remote  section.  Such  a  decrease  of  the  relative  stagna¬ 
tion  temperature  cannot  produce  a  great  influence  on  heat 

transfer . 

6.  By  interaction  of  the  vortex  with  the  spherical  blunt  bo¬ 
dy  a  vortex  breakdown  occurs.  The  vortex  breakdown  is  follo¬ 
wed  by  its  expansion,  formation  of  a  dead  zone  in  front  of 
the  blunt  body  and  flow  turbulation.  This  flow  is  similar  to 
the  flow  in  front  of  a  spiked  body. 

7.  The  vortex  which  is  coaxial  with  the  spherical  blunt  body 
causes  a  heat  flux  decrease  in  the  vicinity  of  the  symmetry 
axis.  At  the  end  of  the  dead  zone  in  the  vicinity  of  the  at¬ 
tachment  line  the  heat  flux  increases  almost  doubles  as  com¬ 
pared  with  the  heat  flux  in  the  stagnation  point  of  the 
sphere  in  the  absence  of  the  vortex.  This  effect  should  be 
taken  into  account  in  analysis  of  heating  bodies  interac 

ting  with  -the  vortex. 

8.  Taking  into  account  the  drastic  influence  of  the  vortex 
on  the  heat  transfer  it  is  reasonable  to  investigate  the  in- 
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ction  of  the  vortex  with  a  cylinder  and  a  plate  at  dif 

^  4.v,^  i/,-.y-i-ov  ?ivis  and  the  streamlined 

^gj-ent  angles  between  the  vortex  axis  anu 

body  surface. 

The  theoretical  model  describing  viscous  diffusion  of  the 
^isymmetric  trailing  vortex  in  compressible  flow  was  sug¬ 
gested  and  developed.  The  main  problem  of  the  model  concerns 
orrelation  of  the  profiles  in  the  initial  section.  The  pro¬ 
files  of  axial  velocity,  temperature  and  circulation  were 
chosen  on  the  basis  of  the  classical  model  of  the  vortex  in 
compressible  gas:  the  flow  was  considered  as  isoenergetic 

and  isoentropic  with  constant  circulation.  At  these  conditi¬ 
ons  the  vortex  contains  vacuum  core  in  the  initial  section. 
The  experimental  results  demonstrate  that  the  vortices  in 
supersonic  flow  really  contain  strongly  decompressed  cores. 

10.  under  action  of  viscous  diffusion  the  vacuum  core  is  fil¬ 
led  by  gas  downstream  of  the  initial  section.  This  process 
was  thoroughly  analyzed  and  self-similar  solutions  near  the 
boundary  of  the  vacuum  core  were  obtained.  Two  types  of  so¬ 
lutions  were  found:  the  solution  with  constant  axial  veloci¬ 
ty  and  another  one  with  constant  total  enthalpy  in  the  ini¬ 
tial  section.  Numerical  calculations  demonstrated  strong 

.influence  of  the  layer  on  the  flow  inside  the  vortex. 

11.  The  asymptotic  expansions  of  the  solutions  far  downstre¬ 
am  from  the  initial  section  were  also  found.  These  solutions 
are  a  generalization  for  compressible  gas  of  well  known  so¬ 
lution  obtained  by  Batchelor  (1964)  for  the  vortex  in  in¬ 
compressible  fluid.  The  asymptotic  and  numerical  results  we¬ 
re  compared  and  good  agreement  between  them  was  demonstra¬ 
ted. 

12.  Systematic  numerical  calculations  of  the  vortices  regro¬ 
uping  under  action  of  viscosity  in  supersonic  flow  were  per¬ 
formed  on  the  basis  of  the  suggested  theoretical  model.  Qua- 


•tatlve  agreement  of  the  numerical  and  experimental  results 
litativ  g  of  the  total  pressure  and 

demonstrated  for  profiles  of 

stagnation  temperature. 

rn-f  the  total  pressure  between  initial  and 
13.  Strong  change  of  experiment  cannot  be  exp- 

r:::":rartrorof  morerular  vlscouse  forces.  -perimental 
Tta  indicate  very  high  level  of  diffusion  in  radial 

tion. 

14  TWO  factors  might  strongly  increase  the  diffusion  of  the 

discontinuity  in  the  vortex  which  was  not  taken 
within  axisymmetrical  model  considered. 

15  Taking  into  account  the  discrepancy  between  theoretical 

■  ntal  results  it  is  worthfull  to  perform  additio- 
and  experimental  results  xi-  int-iude 

.  .  citi.dv  of  the  vortex  structure  and  include 

nai  exoerimental  study  oi  tue 

melsureLnts  of  fluctuations.  This  study  will  °  ^  ; 

mate  turbulent  diffusion  and  develop  further  the 
cal  model  of  the  vortex  in  supersonic  flow. 
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Appendix 


Tables  of  the  measurements  results 


Table  1 


Run  N37  ^  I _ ^ ^ _ L 

Aerator  N1  Alpha  =  2^0°,  Ps/Po  =  f(Phi7)^) 

X  =  2mm _ _ 


qr. 

0 

20 

40 

60 

72 

76 

80 

-45 

0.232 

0.226 

0.228 

0.237 

0.314 

0.305 

0.285 

-30 

0.263 

0.26 

0.261 

0.272 

0.355 

0.338 

0.318 

-15 

0.271 

*  0.268 

*  0.27 

0.28 

0.37 

0.352 

0.325 

0 

*  0.272 

0.267 

0.269 

*  0.281 

*  0.372 

*  0.355 

*  0.328 

15 

0.271 

0.266 

0.268 

0.28 

0.369 

0.35 

0.323 

30 

0.253 

0.241 

0.245 

0.255 

0.339 

0.32 

0.291 

45 

0.191 

0.18 

0.177 

0.187 

0.248 

0.229 

0.202 

90 

95 

96 

97 

98 

99 

100 

-45 

0.163 

0.087 

0.06 

*  0.034 

0.02 

0.016 

0.015 

-30 

0.173 

0.088 

*  0.06 

0.032 

0.02 

0.02 

0.017 

-15 

*  0.179 

*  0.088 

0.057 

0.029 

0.02 

0.019 

0.017 

0 

0.177 

0.085 

0.052 

0.028 

*  0.021 

*  0.019 

*  0.017 

15 

0.174 

0.077 

0.045 

0.022 

0.02i1 

0.018 

0.016 

30 

0.147 

0.057 

0.031 

0.019 

0.02 

0.017 

0.015 

45 

0.093 

0.034^ 

0.019 

0.014 

0.016 

o.oi^ 

0.013 

102 

103 

104 

105 

106 

107 

-45 

0.01 

0.01 

0.012 

0.017 

0.024 

0.031 

-30 

0.013 

0.013 

0.017 

0.027 

0.042 

0.054 

-15 

0.016 

0.019 

0.03 

0.051 

0.075 

0.091 

0 

0.019 

0.026 

0.045 

0.072 

0.095 

0.118 

15 

0.019 

0.029 

0.05 

0.075 

0.102 

*  0.124 

30 

0.021 

0.032 

*  0.054 

0.079 

0.099 

0.122 

45 

*  0.022 

*  0.033 

0.053 

*  0.081 

*  0.107 

0.123 

110 

117 

120 

-45 

0.039 

0.03 

0.031 

-30 

0.068 

0.051 

0.055 

-15 

0.109 

0.079 

0.082 

0 

0.137 

0.091 

0.091 

15 

*  0.142 

*  0.093 

*  0.093 

30 

0.138 

0.093 

0.092 

45 

0.138 

:  0.091 

0.093 

Table  5 

Sphere 

Run  28 

M=3.0  Re=6. 70*10^ 

Tq  =  413  K 

Po=4. 834*10^  Pa  T=1.0  s 

Thermo¬ 

couple 


1.508 
1 . 340 
1 . 173 
1.005 
0.838 
0 . 670 
0.503 
0 . 335 
0 . 168 
0.000 
0 . 084 
0.168 
0 . 251 
0 . 335 
0 .419 
0 . 503 
0.586 
0 . 670 
0.754 
0 . 838 
0 . 922 
1 . 005 
1.089 
1 .173 

1.257 
1 . 340 

1 . 424 
1 . 508 
1 . 591 
1.758 
1 . 925 
2.091 

2 . 258 

2.425 
0 . 168 
0 .335 
0 . 503 
0 . 670 
0 . 838 
1 . 005 
1 .173 
1 . 340 
1.508 


301.7 

306 . 3 

309 . 4 

313.7 

319.5 

325.6 

331.4 
336 . 0 
338.2 

340.8 

339 . 8 

337.6 

333.9 
333.0 

332 . 1 

327.7 

326.2 

322 . 1 

317.6 
315.0 

312.7 

307.8 

308.2 
306 . 0 

302.4 

300.5 

299 . 7 

298 . 6 
299 . 0 

297.8 
298 . 1 
299 . 0 

299 . 9 
301 . 0 
338.4 

330.7 

330 . 9 

325.9 

314 . 6 

309 . 7 

303 . 9 

297 . 9 
292.0 


"^sm  0  ^ 

kW/m"^ 

301 . 7 

7 . 595 

306 . 2 

14.44 

309 . 3 

21.54 

313.7 

29 . 32 

319 . 5 

37 . 19 

325 . 6 

45.94 

331 . 3 

50 . 09 

336 . 0 

50 . 83 

338 . 3 

55 . 36 

340 . 9 

52.59 

339 . 9 

56 . 29 

337 . 7 

56.86 

334 . 0 

57 . 60 

333 . 0 

55.81 

332 . 2 

51 . 80 

327 . 8 

52 . 32 

326 . 2 

49 . 82 

322 . 1 

50.59 

317 . 6 

44 . 19 

315 . 0 

38.90 

312 . 8 

35 . 17 

308 . 0 

30 . 41 

308 . 3 

26 . 29 

306 . 0 

22 . 86 

302 . 3 

17 . 87 

300 . 6 

13 . 84 

299 . 6 

11 . 32 

298 . 6 

8 . 158 

299 . 0 

5.534 

297 . 8 

4.710 

298 . 1 

4 .404 

299 . 0 

4 . 139 

301 . 0 

4 . 128 

338 . 5 

54 .14 

330 . 8 

56 . 78 

330 . 8 

50.72 

325 . 9 

46 . 34 

314 . 6 

41 . 28 

309 . 7 

32.77 

303 . 9 

25.99 

297 . 8 

16.09 

292 . 0 

11 . 33 

76 . 95 

73 .79 

71 . 63 

68 .63 

64 . 62 
60 . 39 

56 .41 
53 . 19 

51 . 57 

49 . 76 

50 . 44 
52 . 02 

54 . 58 
55 . 22 
55 . 82 
58 . 85 
59 . 93 

62 . 78 
65 . 88 
67 . 69 
69 . 26 

72 . 59 

72 . 36 

73 . 96 
76 . 48 
77 .71 

78 . 37 
79 . 04 

78 . 80 

79 . 63 

79 . 42 

78 . 78 
78 . 17 

77 . 42 

51 . 44 

56.77 
56 .73 
60 . 16 
68 . 02 

71 . 38 

75 . 42 

79 . 63 

83 . 63 


0.0987 
0 . 1957 
0 . 3007 
0 . 4272 
0.5755 
0.7607 
0.8879 
0.9556 
1 . 0735 
1.0569 
1 .1160 
1.0932 
1 . 0552 
1.0107 
0 . 9279 
0 . 8891 
0.8312 
0 . 8058 
0.6707 
0.5748 
0 . 5077 
0.4190 
0.3633 
0 . 3091 
0.2337 
0 . 1781 
0 .1445 
0 . 1032 
0 . 0702 
0 . 0592 
0 . 0555 
0 . 0525 
0.0528 
0 . 0533 
1 . 0524 
1.0002 
0 . 8940 
0 . 7703 
0 . 6069 
0 . 4590 
0 . 3446 
0 . 2020 
0 .1354 


sphere  with  vortex  generator  O  20"  Run  2V 

,,3.0  Re=6. 86.10=  T„  =  410  K  P„  =4,889.10=  Pa  1=1.0  s 


Table  6 
Run  29 


vi=3 .0 

Re=6. 86*10 

Xhermo- 

couple 

S/R 

1 

1 . 508 

2 

1 . 340 

3 

1.173 

4 

1 . 005 

5 

0 . 838 

6 

0.670 

7 

0 . 503 

8 

0.335 

'^sm  0  ^  I 


9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


0.168 
0 . 000 
0 . 084 
0 . 168 
0 . 251 
0 . 335 
0 . 419 
0 . 503 
0.586 
0 . 670 
0.754  I 
0 . 838 
0 . 922 
1 . 005 
1 . 089 
1.173 

1.257 
1 . 340 

1.424 
I  1.508 

1 . 591 
1.758 
1.925 
2 . 091 

2.258 

2 .425 
0.168 
0 .335 
0.503 
0.670 
0.838 
1.005 
1 . 173 
1  .  340 
1  p;nR 


311 . 8 

320 . 3 

329.7 

339 . 3 

347 . 3 

354 . 3 

351 . 3 

338.4 

331.9 

335 . 1 

335 . 2 

335.4 

335 . 8 

339 . 6 

346 . 1 

345 . 9 

346 . 2 

343 . 2 

339.4 

336 . 4 

332.7 

328 . 9 

326.8 
j  321.8 

317.2 

313.5 

310.5 
308.0 

306 . 8 

306 . 1 

306 . 9 

308 . 9 

310.9 

312.6 

332.4 
326.0 

324 . 3 
321 . 0 
314 . 0 

310 . 3 

305 . 3 

301.6 

297.2 


kW/m^  kW/m 


311.8 

320 . 3 

329.7 

339 . 4 

347.4 

354 . 3 

350.9 

338.2 

331 . 9 
335.0 
335 . 0 

335.5 
336 . 0 

340.1 

346 . 5 

346 . 3 

346.6 

343.5 

339.8 

336.7 
333.0 

329.1 
j  326.9 

322 . 0 

317 . 2 

313.6 
310.5 
308 . 0 

306 . 8 

306 . 1 
307.0 
309 . 0 
311.0 

312.7 
332.5 
326.0 

324 . 3 

321 . 1 

314 . 3 

310 . 4 

305 . 5 

301 . 7 
297.4 


11 .94 
23 . 21 

34.13 
43.72 
52.62 

66 . 13 
72 . 61 
56 . 54 

43.67 
39 . 06 

42.97 
48 . 31 
59.03 
68 . 64 
70 . 24 
79 . 17 

77.46 

80.06  I 

69 . 69 

61.13 
53.48 

45.95 
38.77 

33 . 51 
25.73 
19 . 35 

15 .46 
10 .42 
6 . 367 
5 . 322 
6.207 
7 .951 
9.814 
11.03 
37 . 61 
41.17 

41.68 
44 . 99 
51 . 10 

47 . 29 
39 . 27 

24 .13 

17 . 13 


68.07 

62 . 14 
55 . 60 
48.82 

43 . 30 

38 .51 
40 . 84 

49 . 68 
54 . 04 
51 . 88 
51.92 
51.54 
51 . 21 
48 . 39 
43 . 91 
44.09 
43 . 84 
45.99 
48 . 59 

50 . 69 
I  53.29 

55 . 98 
57 . 50 

60 . 95 

64 . 30 
66 . 80 
68 . 92 

70 . 69 

71.52 
72 . 00 
71.40 
70 . 02 
68 . 63 
67 .44 
53.65 
58.16 
59 . 31 
61.59 
66 . 33 

68 . 99 

72.46 
75 . 05 
78 . 09 


q/qo 


0.1755 
0.3736 
0.6139 
0 . 8956 
1 . 2152 
1.7171 
1.7780 
1 .1382 
0 . 8081 
0 . 7528 
0 . 8277 
0 . 9373 
1 . 1527 
1 . 4186 
1 . 5997 
1 .7956 
1 .7669 
1.7406 
1 . 4344 
1 . 2059 
j  1.0035 
0 . 8209 
0 . 6742 
0 . 5498 
0 . 4002 
0 . 2897 
0 . 2243 
0 . 1474 
0 . 0890 
0.0739 
0  -  0869 
0.1136 
0 . 1430 
0 . 1635 
0.7010 
0 .7079 
0 .7027 
0 .7305 
0.7704 
0.6855 
0 . 5420 
0.3215 
0 . 2193 


Table  8 
Run  43 


‘sphere  with  vortex  generator  a  30  Run 

,.3.0  Re.6. 76*10=  To*  414  K  P„=4. 891*10=  Pa  1-1.0  s 


M=3.0 

Re=6. 76*10 

Thermo¬ 

couple 

S/R 

1 

1 . 508 

2 

1 . 340 

3 

1.173 

4 

1.005 

5 

0.838 

6 

0.670 

I  7 

0 . 503 

8 

0 . 335 

0 . 168 
0 . 000 
0 . 084 
0 .168 
0 . 251 
0.335 
0.419 
0 . 503 
0 . 586 
0 . 670 
0.754 
0 .838 
0.922 
1 .005 
1 .089 
1 . 173 

1 . 257 
1 . 340 

1 .424 
1.508 
1 .591 
1.758 
1.925 
2.091 

2.258 

2.425 
0.168 
0.335 
0.503 
0.670 
0 . 838 
1 . 005 
1 . 173 
1 . 340 
1 . 508 


299 . 9 

304 . 5 
309 . 0 

312.9 

317.6 

322 . 2 

324 . 4 

329 . 3 

327 . 4 

330 . 1 

328 . 4 

327.7 

326.7 

323 . 2 

320 . 3 

318.7 

316.8 

314 . 3 

306 . 9 

307 . 2 

306 . 3 

301 . 7 

303 . 8 

300 . 6 

298 . 7 

297 . 9 

299 . 9 

296.5 

298 . 1 

299 . 6 

302 . 4 
j  303.5 

306 . 6 
306.0 
329 . 0 

335 . 2 

333.2 

318.3 

304 . 5 

305 . 8 

309 . 8 

308 . 5 

300.4 


Tsmo  ^ 


299 . 8 

304 . 5 
309 . 0 

312.8 

317.5 

322 . 1 

324 . 1 

329 . 3 

327 . 2 
330.0 

328 . 4 

327.6 

326 . 5 

323 . 2 

320.3 

318.8 

316.7 

314 . 3 

306 . 8 

307 . 2 

306 . 2 

301 .7 

303.7 

300 . 5 

298.6 

297.8 

299 . 8 

296.4 
298 . 0 

299 . 5 

302.2 

303 . 3 

306 . 4 

305 . 8 

328 . 3 

333 . 5 
333 . 0 

318.5 

304 . 7 

305 . 5 

309 . 5 

308 . 4 

300 . 4 


q  ^  „  qo  '  2 

kW/m^  kW/m'" 


6.333 

11.72 

16 . 97 
22 . 22 

27 . 82 
33 .70 

35.98 
35.91 
38.27 

35.72 
36 . 88 

35 . 38 

35 . 39 

35 . 26 

33 . 26 

34 . 65 

33.55 
33 . 97 

29.65 

26 . 17 

23 . 36 
20 . 11 
17 . 33 
15.05 
12 . 04 

9 . 569 
7 . 845 
5 . 914 
4 . 154 
3.943 
4 .645 

5 . 570 
6 . 679 
7 . 315 
34.49 

52 .51 

52 . 36 

39.51 
24.78 

22 . 82 
27.75 
20.64 
16.59 


79 . 54 
76 . 30 
73 . 16 

70 . 51 

67.21 
64.05 

62 . 61 

58.99 
60 . 44 

58.55 

59 . 62 

60 . 21 
60 . 93 
63 . 28 
65.25 
66 . 35 
67.75 
69 . 48 
74.67 
74.41 
75 . 10 
78 . 24 

76 .82 
79 . 08 
80 . 38 

j  81.00 

79 . 55 
81 .93 
80 .81 
79.79 
77 . 90 
77 .11 
74 . 98 

15. 37 

59.72 
56 . 07 
56 .43 
66 . 50 

76.18 

75 . 62 
12. IZ 
73 . 58 

79.18 


q/qo 


0 . 0796 
0 . 1536 
0 . 2320 
0 . 3152 
0.4139 
0 . 5262 
0 . 5747 
0 . 6087 
0.6332 
0 . 6102 
0.6185 
0 . 5875 
0 . 5809 
0 . 5571 
0 . 5097 
0.5222 
0.4952 
0 .4889 
0 . 3971 
0 . 3516 
0 . 3110 
0.2570 
0 . 2256 
0 . 1903 
0.1498 
1  0.1181 
0 . 0986 
0 . 0722 
0 . 0514 
0 . 0494 
0 . 0596 
0 . 0722 
0 .0891 
0.0971 
0 . 5774 
0 .9365 
0 . 9279 
0 . 5942 
0 . 3253 
0 .3018 
0 . 3813 
0 . 2805 
0.2095 


0.90 


Generator  N2 


Test  M24 


X=100  mm 


Run  N37:  X=2mm 


-20  0  20  40  60 

Phi,  degrees 


Run  N36:  X=110mm 


-20  0  20  40  60 

Phi,  degrees 


Run  N38:  X=2mm 


■ —  Y'=40mm 
□ —  Y’  =  100nnm 
—  Y'=106mnn 


Run  N41:  X=110mm 


Phi,  degrees 


to,  degrees  C 


x=2mm 


to,  degrees  C 


x=llOmm 


Stagnation  temperature 
after  generator  N1 , 
a=  20  ° 


Fia.  31 


theory 


Sphere  with  generator 


h(o. 2f  Xg/D=l - 17 


Run  30 


theory 
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